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Abstract 



The present theoretical knowledge about photons and hard photoproduction pro- 
cesses, i.e. the production of jets, light and heavy hadrons, quarkonia, and prompt 
photons in photon-photon and photon-hadron collisions, is reviewed. Virtual and 
polarized photons and prompt photon production in hadron collisions arc also dis- 
cussed. The most important leading and next-to-leading order QCD results are 
compiled in analytic form. A large variety of numerical predictions is compared to 
data from TRISTAN, LEP, and HERA and extended to future electron and muon 
colliders. The sources of all relevant results are collected in a rich bibliography. 
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I. INTRODUCTION 



The photon is a fascinating particle. It is the fundamental gauge boson of Quantum 
Electrodynamics and as such one of the best studied elementary particles of the Standard 
Model of particle physics. However, in high-energy reactions the photon exhibits a complex 
hadronic structure, which is far less well understood. Investigations took off with fixed-target 
experiments and soft vector-meson production more than 40 years ago, which established 
the Vector Meson Dominance model of fluctuations between photons and vector mesons. 
With the advent of Quantum Chromodynamics (QCD) as the theory of strong interactions 
in the 1970s, the pointlike nature of the photon, its coupling to quarks, and the perturbative 
evolution of the photon structure moved to the center of interest. For some time it seemed 
possible to actually calculate the hadronic structure of the photon. Unfortunately it turned 
out shortly later, that this was only possible at asymptotically large energies, which could 
not be reached with the available accelerators. In the mid-80s, the photon beam energies 
were high enough to prove that there were two different types of photon interactions. The 
photon could interact directly with the quarks and gluons in the hadronic target, but it 
could also resolve into a hadronic structure and the partonic constituents of the photon 
could participate in the hard scattering. At about the same time, the two processes were 
found to be theoretically related. The pointlike photon structure was shown to be singular 
and in need of regularization by a non-perturbative hadronic boundary condition, which had 
to be determined from theoretical models or deep-inelastic electron-photon experiments at 
the e+e" colliders PETRA and PEP. 

During the last 15 years, a tremendous experimental and theoretical effort has been 
invested into refining the picture of the photon, testing its structure, and using the pho- 
ton as a tool for studying the production and properties of jets, light and heavy hadrons, 
quarkonia, and prompt photons. Experimentally this was facilitated by the construction 
and operation of the e + e~ colliders TRISTAN and LEP with center-of-mass energies be- 
tween 58 and 210 GeV and the HERA ep collider with a center-of-mass energy of 300 GeV, 
which has just been upgraded for high-luminosity operation. The experimental results on 
the structure and interactions of the photon have been discussed at many general and the 
topical "PHOTON" conferences, and they have recently been reviewed by several authors 
(Erdmann, 1997; Abramowicz and Caldwell, 1999; Nisius, 2000; Krawczyk, Staszel, and 
Zembrzuski, 2001). The experimental advances were paralleled by similar theoretical im- 
provements. Until the mid-80s, determinations of the photon structure and calculations of 
photoproduction processes were performed only in leading order (LO) of perturbative QCD, 
and a review of the LO formalism exists for prompt photon, jet, and particle production in 
hadron collisions (Owens, 1987). Since then, several photon parameterizations and many 
photoproduction calculations have been performed in next-to-leading order (NLO). In com- 
parisons with experimental data, these calculations have provided stringent tests on QCD 
and precise information about the photon structure and fragmentation and about the forma- 
tion of jets and hadrons (for short reviews see Klasen, 1997a, 1999b; Kniehl, 1997; Kramer, 
1996, 1998, 1999). 

While LO predictions are straightforward to calculate, they are unfortunately not very 
precise and can only be used to estimate hard photoproduction cross sections within factors 
of two. The reason for this is that LO QCD cross sections depend strongly on the renormal- 
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ization scale in the strong coupling and on the factorization scales in the photon and proton 
parton densities, which are commonly varied by factors of two around the physical scale 
to estimate the relative theoretical error. Specific problems exist for jet and heavy-flavor 
production. In LO every jet corresponds to just one parton, and jet algorithms cannot be 
implemented. Only beyond LO this becomes possible by combining two or more partons into 
a single jet. Heavy flavors can be treated as massive or massless particles, and only in NLO 
the logarithmic mass and scale dependences in fixed and variable flavor number schemes can 
be made explicit. The large luminosities and trigger rates at modern high-energy colliders 
and detectors lead to great improvements in the statistical accuracy of the experimental 
data, which exceeds the theoretical precision of LO calculations. In NLO the uncertainties 
can be reduced to a reasonable level (a few percent). This precision is generally limited to 
scattering processes with four external legs, and it has now been reached for almost all pho- 
toproduction processes. Exceptions are the real photoproduction of three jets, two hadrons 
or photons, and various quarkonium states, and the virtual or polarized photoproduction of 
hadrons and prompt photons. In some cases NLO precision may be insufficient to describe 
the experimental data as will be discussed later. It would be necessary to work at next-to- 
next-to-leading order, but this precision has so far been feasible only for inclusive processes 
like the total rates of e + e~ and ep scattering into hadrons. 

It is the aim of this Review to give a complete theoretical description of hard photopro- 
duction processes, ranging from the generation of real, slightly virtual, and also polarized 
photons with lepton beams, over the current knowledge about the photon structure, to the 
methods and applications of NLO QCD calculations for jet, hadron, and prompt photon 
photoproduction. While the theoretical tools cannot be described in full detail, the main 
lines of argument and techniques are explained conceptually, and the most important re- 
sults are stated explicitly. The predictions are compared only to the most recent and precise 
data from TRISTAN, LEP, and HERA in order to reach up-to-date conclusions about the 
applicability of QCD and the properties of jets, hadrons, and photons. Selected predictions 
are made for future accelerators like linear e + e~ and circular muon colliders or a future ep 
machine in order to provide a look ahead, and a rich bibliography guides the reader to all of 
the relevant theoretical and experimental literature. In this sense, this Review will hopefully 
serve as a compendium of the current state of the art in hard photoproduction. 
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II. PHOTON SPECTRA 



In the first generation of photoproduction experiments, real photons with energies below 
60 GeV were generated by pion decays or electron bremsstrahlung and scattered off nu- 
clear targets in order to study soft particle production and total cross sections (Paul, 1992). 
Measurements of hard photoproduction of photons and mesons began with the CERN ex- 
periment NA14 at beam energies between 50 and 150 GeV (Auge et al., 1986). Later 
energies up to 400 GeV were reached at Fermilab. They led to a first observation of jets 
in photoproduction (Adams et al, 1994), but were still too small for definite tests of QCD. 
Center-of-mass energies of 200 to 300 GeV have recently been reached at the e + e~ colliders 
PETRA, PEP, TRISTAN, and LEP and at the ep collider HERA. Here spacelike, almost 
real bremsstrahlung photons are exchanged during the hard collision. At future linear e + e~ 
colliders large particle bunch densities are needed to reach high luminosities. Then addi- 
tional beamstrahlung photons will be created before the hard interaction by the coherent 
action of the electromagnetic field of one bunch on the opposite one. If the electron beams 
are collided with additional high-energy laser beams, real photons can be produced through 
Compton scattering. Thus three different mechanisms can contribute to photon scattering 
at high-energy colliders: bremsstrahlung, beamstrahlung, and laser backscattering. In this 
Section, the corresponding photon energy spectra will be discussed. 

A. Bremsstrahlung 

If the outgoing and incoming leptons in a hard scattering process are almost collinear, 
the calculation of the corresponding cross section can be considerably simplified by using 
the Weizsacker- Williams or Equivalent Photon approximation (for a review see Budnev et 
al, 1974). Current conservation and the small photon virtuality lead to a factorization of 
the lepton scattering cross section into a broad-band photon spectrum in the lepton and 
the hard photon scattering cross section. Already in the 1920s, Fermi (1924) discovered the 
equivalence between the perturbation of distant atoms by the field of charged particles flying 
by and the one due to incident electromagnetic radiation. His semi-classical treatment was 
then extended to high-energy electrodynamics by Weizsacker (1934) and Williams (1934) 
independently, who used a Fourier analysis to unravel the predominance of transverse over 
longitudinal photons radiated from a relativistic charged particle. In the fifties, Curtis 
(1956) and Dalitz and Yennie (1957) gave the first field-theoretical derivations and applied 
the approximation to meson production in electron-nucleon collisions. 

In the Weizsacker- Williams approximation, the energy spectrum of the exchanged pho- 
tons is given by 



The subleading non-logarithmic terms (Kessler, 1975; Frixione et al, 1993) modify the 
cross section typically by 5%. a = 1/137 is the electromagnetic fine-structure constant, 




a \l + (l-x)\ QLx(l-z) 



2n x mfx 2 




(2.1) 
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x = E 1 /Ei (I = e, fi) is the energy fraction transferred from the lepton to the photon, mi 
is the lepton mass, and Q^ax = E? (1 — x ) * s the maximal photon virtuality for lepton 
scattering angles below 9 ma _ x . This angle can be determined by tagging the outgoing lepton in 
the forward direction or by requiring that it is lost in the beam pipe (anti-tagging). When no 
information about the scattered lepton is available, one has to integrate over the whole phase 
space, thus allowing large transverse momenta and endangering the factorization property of 
the cross section. Bremsstrahlung photons can be generated by electron and muon beams, 
but for the latter the photon density is smaller by approximately a factor of two due to the 
larger muon mass (Klasen, 1997b). 



B. Beamstrahlung 



Future circular e + e colliders with center-of-mass energies above y/S = 500 GeV would 
suffer from very high synchrotron radiation. They must therefore have a linear design and 
dense particle bunches in order to still obtain large luminosities. Inside the opposite bunch, 
electrons and positrons experience transverse acceleration and radiate beamstrahlung. The 
spectrum is controlled by the beamstrahlung parameter 



T 



5r 2 e E e N 



Qaa z (a x + o y )m e ' 



(2.2) 



which is proportional to the effective electromagnetic field of the bunches and depends on 
the classical electron radius r e = a/m e = 2.818 • 10~ 15 m, the beam energy E e , the total 
number of particles in a bunch N, and on the r.m.s. sizes of the Gaussian beam a x , a y , a z . 
For not too large Y (Y < 5), Chen (1992) derived the approximate spectrum 



/ 7 b ;r(*) = 




X 3(1 — x) s e 



-2x/[3T(l-a;)] 



X 



+ \hrr 



1 - 



g(x)N, 



with 



9(z) = 1 - \ 



1 - :r)/l - I" + i 



x 



(1-X)3. 



The average number of photons radiated per electron throughout the collision is 

5a 2 a 7 m P T 



-r, /•:, \ 1 -I T" 



(2.3) 



(2.4) 



(2.5) 



Current design parameters for possible future e + e colliders are listed in Tab. I, and the 
spectra corresponding to the v^S* = 500 GeV designs are shown in Fig. I. This figure also 
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Collider 


TESLA 


ILC-A 


ILC-B 


ILC-C 


CLIC 


Last update 


3/01 


4/00 


4/00 


4/00 


3/01 


VS (GeV) 


500 


535 


515 


500 


500 


Particles/Bunch (10 10 ) 


2 


0.75 


0.95 


1.1 


0.4 


a x (nm) 


553 


277 


330 


365 


202 


o v (nm) 


5 


3.4 


4.9 


7.6 


2.5 


a z (jum) 


300 


90 


120 


145 


30 


T 


0.05 


0.14 


0.11 


0.09 


0.28 


V<5 (GeV) 


800 


1046 


1008 


978 


1000 


Particles/Bunch (10 10 ) 


1.4 


0.75 


0.95 


1.1 


0.4 


cr^ (nm) 


391 


197 


235 


260 


115 


a y (nm) 


2.8 


2.7 


3.9 


5.4 


1.75 


cr z (//m) 


300 


90 


120 


145 


30 


T 


0.10 


0.38 


0.30 


0.25 


0.98 



TABLE I. Current design parameters for possible future linear e + e colliders (Brinkmann et al., 
2001; Adolphsen et al, 2000; Kamitani and Rinolfi, 2001). 



displays the bremsstrahlung spectra for electrons and muons, which have been integrated 
over the photon virtuality up to an upper bound Q 2 nax = AE?(1 — x) for untagged outgoing 
leptons. Both the bremsstrahlung and beamstrahlung spectra are soft, but for particular 
collider designs beamstrahlung can be more important than bremsstrahlung over a wide 
range in x. Beamstrahlung is completely negligible at muon colliders, since the parameter 
T M = m^/m^Te, which controls the spectrum mostly through the first exponential function 
in Eq. (2.3), is about five orders of magnitude smaller than for electrons (Klasen, 1997b). 



C. Laser backscattering 

Real photons of very high energy can be produced, if laser photons are backscattered off 
electrons. The laser backscattering spectrum (Ginzburg et al, 1984) 



where 



and 



+ 



N c + 2\ e P c N' c 



l-x + 
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l-x X(l-x) 
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Photon Spectra at Lepton Colliders with Vs = 500 GeV 
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FIG. 1. Photon energy spectra at lepton colliders with = 500 GeV. The bremsstrahlung 
spectrum depends logarithmically on \/~S and is roughly two times larger for electrons than for 
muons. The electron beamstrahlung spectra depend strongly on the design parameters and rise 
with y/~S. The laser backscattering spectrum is completely independent of V~S, if the the center-of- 
mass energy of the electron-laser photon collision is kept fixed, but it depends on the polarizations 
of the electron A e and of the laser photon P c . 



N' 



c 



1 + ]n fi + X )-- + — — L 

XJ V ; 2 1 + X 2 1 + 



xy 



(2.8) 



are related to the total Compton cross section, depends on the center-of-mass energy of the 
electron-laser photon collision s ei through the parameter X = s ei jvr? e — \. The optimal value 
of X is determined by the threshold for the process 77 — > e + e~ and is X = 2(1 + y/2) ~ 
4.83 (Telnov, 1990). If this value is kept fixed, the laser backscattering spectrum becomes 
independent of \/~S. A large fraction of the photons is then produced close to the kinematic 
limit x < x max = Xj (1+X) ~ 0.828, so that one obtains an almost monochromatic "photon 
collider." The monochromaticity of the produced photons can be improved further if the 
helicities of the electron A e and of the laser photon P c satisfy the condition 2\ e P c = — 1 
(see Fig. 1). For electron beams, the optimal pulse energies (~ 1 J) and repetition rates 
(~ 1 kHz) can be provided by current laser technology. For muon beams, one would need 
much higher flash energies (~ 1 GJ), which are currently not feasible (Klasen, 1997b). High- 
energy interactions of laser-backscattered photons have yet to be experimentally observed, 
and the simulation tools for realistic photon spectra including low-energy tails and non-linear 
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effects have to be refined in order to reliably estimate the luminosities and event rates for 
photon colliders. 
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III. PHOTON STRUCTURE 

The first generation of photoproduction experiments on nuclear targets revealed striking 
similarities with purely hadronic processes like pion scattering. The total cross sections 
showed strong resonances below center-of-mass energies of 3 GeV and flat behavior above, 
and the elastic scattering cross sections fell strongly with the momentum transfer (Bauer et 
ai, 1978). The most striking feature, however, was the copious production of vector mesons, 
particularly of the p meson, which led Stodolsky (1964) to postulate an analogy between the 
isovector electromagnetic current and the p field operator and a relation between diffractive 
7P and elastic pp cross sections (Ross and Stodolsky, 1966). This analogy was subsequently 
extended in the Vector Meson Dominance (VMD) model, where the photon was viewed as 
a superposition of p, u, and <f) mesons with a small, if not negligible, pointlike contribution 
(Stodolsky, 1967; Joos, 1967). In the Generalized Vector Dominance (GVD) model (Sakurai 
and Schildknecht, 1972), all mesons, which carried the same quantum numbers as the photon 
(J PC = 1 ), were included as contributions to the photon structure, e.g. also the heavier 
J/\l/ meson. 

The hadronic picture of the photon was revolutionized with the advent of QCD. Witten 
(1977) showed, that the photon structure had an anomalous (pointlike) component, which 
could be understood as a short-time fluctuation into quark- ant iquark pairs and gluons and 
could be calculated perturbatively. QCD also predicted direct photon scattering off quarks 
and gluons, which was subsequently observed in fixed-target experiments at higher ener- 
gies and momentum transfers (Auge et ai, 1986). Thus it became clear that the photon 
had a pointlike and a hadronic component, that they were important at different momen- 
tum transfers, and that they induced photon and parton scattering reactions, respectively 
(Brodsky et ai, 1978, 1979; Lewellyn-Smith, 1978). Much of the present knowledge about 
the hadronic structure of the photon has been obtained from measurements of the photon 
structure functions and Fl in the deep-inelastic scattering process e{k)^{p) — > e(k')X 
at e + e~ colliders with 



where the flux of quasi-real photons is given by the Equivalent Photon approximation, Eq. 
(2.1), and where the kinematic variables Q 2 = —q 2 = — (k — k') 2 , x = Q 2 /(2p ■ q), and 
V — ip ■ o)/(p ' k) are experimentally determined. 

A. Evolution equations 

In NLO of QCD, an initial-state photon can split into a massless quark-antiquark pair, 
which then interacts in the hard scattering, and a collinear singularity is encountered. In di- 
mensional regularization ('t Hooft and Veltman, 1972; Bollini and Giambiagi, 1972a, 1972b; 
Gastmans and Meuldermans, 1973; Marciano and Sirlin, 1975; Marciano, 1975; Leibbrandt, 
1975), where scattering matrix elements and phase-space factors are evaluated in n = 4 — 2e 
dimensions, the singularity manifests itself as a 1/e-pole multiplying the space-like photon- 
quark splitting function P q ^(x) in 




(3.1) 
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T q ^ y (x, Mj) = 5 qi 5(1 — x 



1 a r(l -e) (4nfi 



2\ £ 



e2nT(l-2e) \ M) 
x P q ^(x) + 0(e, a 2 , aa s ) 



5 qi 5(1 — x) — 



a 



x —P q ^(x) + 0(e,a ,aa s ). 
lix 



(3.2) 



x is the longitudinal momentum fraction of the quark in the photon, 7^ = 0.5772 ... is the 
Euler constant, and the scale \x has been introduced to preserve the dimension of physical 
quantities. The factorization theorem (Amati, Petronzio, and Veneziano, 1978a, 1978b; 
Ellis et al, 1979; Collins, Soper, and Sterman, 1988) ensures, that the collinear singularity 
appearing in the transition function T is universal, i.e. independent of the LO scattering 
process, and can be absorbed into a renormalized quark density in the photon 



f qh (x,Mj) = f qh (x)+ T q ^(Mj)®f lh (x) 



(3-3) 



at a factorization scale Mj. Thus logarithmic dependences on the artificial scale Mj are 
induced in / ? / 7 and in the NLO partonic scattering cross section, which cancel up to higher 
orders in the perturbative expansion. In the modified minimal subtraction (MS) scheme 
(Bardeen et al, 1978) no additional finite terms are subtracted, since these depend generally 
on the hard scattering process. In Eq. (3.3) the x-space convolution is defined as 



./ 



1/7 



(*>=/ 



i* 

y 



X 



fih(v), 



(3.4) 



f q /y(x) = 5 qi 5(l — x) + 0(a) is the bare quark density, and f^^(x) = 5(1 — x) + 0(a 2 ) 
is the bare photon density, which gets only renormalized at next-to-next-to-leading order 
(NNLO) in a. A renormalized gluon density in the photon f g / y (x, Mj) is generated in a 
similar fashion by gluon radiation from a quark- ant iquark pair at 0(aa s ), where a s (ji 2 ) 
is the running coupling constant of quarks and gluons. Beyond LO, the direct photon 
scattering processes are thus intimately connected to the resolved processes. In this Section, 
we review the NLO QCD evolution equations, boundary conditions, factorization schemes, 
and hadronic solutions of the photon structure function. 

The evolution of the parton densities in the photon with the scale M 2 can be calculated 
in perturbation theory by taking the logarithmic derivatives of / ? / 7 , f g / 7 , and / 7 / 7 with 
respect to M 2 . This leads to a coupled system of integro-differential equations (De Witt et 
al, 1979) 



d/ g/7 (Q 2 ) 
dlnQ 2 
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2tt q ^ 
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(3.5) 
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P a ^ q ®f qh {Q 2 ) + P^ 9 ®f 9 h(Q 2 ) 



d/ 7/7 (g 2 ) 
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2tt 7< ~ 7 



/ 7 / 7 (g 2 ) + 
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2^ 



x 



P^g ® / 9 / 7 (Q 2 ) + P^g ® fghiQ 2 ) 



These evolution equations differ from the well-known hadronic case by the inhomogeneous 
terms of 0(a), which arise from the pointlike coupling of photons to quarks. Eqs. (3.5) are 
given for a single quark flavor q, but 2Nf light quarks and antiquarks are easily accommo- 
dated by summing over the index q from one to 2Nf whenever it appears twice. The renor- 
malization scale fi and the factorization scale Mf have been identified with the physical scale 
Q. The x-space convolution Pj<-i <E> fi/^iQ 2 ) reduces to a simple product Pj t _ i (n)fi/ 1 (n, Q 2 ) 
in Mellin n-space, where the n th moment is defined as f(n) = Jq dxx n ~ 1 f(x). In moment 
space, analytical solutions to the evolution equations can be found, but they have to be 
transformed back to rr-space for physical cross section predictions. Alternatively, the evo- 
lution equations can be solved directly in x-space by iteration, but this requires a careful 
separation of LO and NLO terms in order to avoid spurious higher order terms. An original 
ansatz by Rossi (1984) was subsequently generalized by Da Luz Vieira and Storrow (1991) 
in order to allow for non-zero input parton densities. 

The LO and NLO splitting functions Pj<-i have been calculated by Curci, Furmanski, and 
Petronzio (1980) for the flavor-nonsinglet case, by Furmanski and Petronzio (1980) for the 
coupled gluon and flavor- singlet cases, and also by Floratos, Kounnas, and Lacaze (1981). 
Here we review only the LO results (Altarelli and Parisi, 1977) 



P q ^ q (x) = C F 
P g ^ q (x) = C F 

P g ^ g (x) = 2N C 
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+ x(l-x)-2 
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—N c - -N f 

6 ° 3 / . 



S(l-x) + 0(a s ), 



(3.6) 



where the -(--distributions are defined as usual with a test function f(x) in the integral 



C dxf(x)g + (x) = f dx[f(x) - f(l)]g(x). 
o Jo 



(3.7) 

The photon-quark splitting function can be obtained in LO from P q <- g by the transformation 
P q ^ = 2N c e 2 q P q ^ g . C F = (N 2 - 1)/(2N C ) and N c = 3 are the SU(3) color factors, 
T R = 1/2, Nf is the number of active flavors, and e q is the fractional quark charge. Since 
there is no direct coupling of photons to gluons, the photon-gluon splitting function enters 
only in NLO. It can be obtained from the NLO gluon-gluon splitting function by replacing 
the appropriate color factors and dropping the part proportional to 6(1 — x) (Fontannaz and 
Pilon, 1992) 
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(3.8) 



B. Boundary conditions and factorization schemes 



The general solutions of the evolution equations Eqs. (3.5) are given by the sums of 
pointlike ("anomalous") and hadronic contributions fi/ 1 {Q 2 ) = ff/^iQ 2 ) + / l y 7 i (Q 2 )- Due to 
the pointlike coupling of the photon to quarks, the former can be calculated perturbatively at 
asymptotically large scales Q 2 and large Bjorken-a; in LO (Witten, 1977) and NLO (Bardeen 
and Buras, 1979). The pointlike solution is of the form 



a 
2k 



a s (Q 2 



-cii + bi + 0(a s ) 



(3.9) 



where Oj and hi are known analytic functions in moment space. Scaling is already violated at 
LO, 0(a/a s (Q 2 )), and the pointlike solution dominates at large Q 2 . As the parton-parton 
cross sections are of 0(a 2 ), the resolved photon contributions are of the same order 0(aa s ) 
as the direct photon-parton contributions. The pointlike solutions turn out to be singular 
at low x and moderate Q 2 (Duke and Owens, 1980) and have to be regularized by boundary 
conditions fffi(Qo) at some low starting scale Qq (Gliick and Reya, 1983). 

By combining the parton distributions functions f q /-y(Q 2 ) and / 9 / 7 (Q 2 ) with the appro- 
priate Wilson coefficients (Bardeen et al, 1978) 



C q (x) = C F 
C g (x) 

and (Bardeen and Buras, 1979) 



x 



x 



In 



1 



x 



X 



+ ^(9 + 5s) 



(3.10) 



T, 



R 



(x 2 + (1 - x) 2 ) In — — + 8x(l -x)-l 
v ' x 



C 1 {x) = 2N c C g {x) 



(x 2 + (1 - x) 2 ) In 



l-x 



X 



+ 8x(l -x) - 1] , 
one obtains the NLO photon structure function in the MS scheme 

r a s (Q 2 



(3.11) 



i?(Q 2 ) = 5> e ; \f qh (Q 2 ) + 



2tt 



+ C g ®f gh {Q 2 ) 



a 



(3.12) 



The factor of two arises from the fact that / g / 7 (<5 2 ) = f q /~/(Q 2 ) due to charge conjugation 
invariance. Gliick, Reya, and Vogt (1992a) observed that the direct term C 7 (x) contains 
a term ln(l — x), which diverges in the large-a; region and is better absorbed in the quark 
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distributions f^{Q 2 ) = f^{Q 2 )+a/{27c)e 2 q C 1 . This also affects the NLO photon splitting 
functions 



pDIS 7 
pDIS 7 



P, 



MS 



e 2 P 
c <? 1 q^q 



a 
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(3.13) 



The definition of the gluon densities in this new DIS 7 factorization scheme remains un- 
changed. Equivalently in the MS scheme one can absorb C 7 (Gordon and Storrow, 1992a, 
1997) or its process-independent part (Aurenche, Fontannaz, and Guillet, 1994a) into point- 
like initial quark distributions. 

Heavy quarks with mass and velocity squared /3 2 = 1 — 4m\xj [(1 — x)Q 2 ] contribute 
to the photon structure function through the Bethe-Heitler process Y(Q 2 )l ~~ hh with 
(Budnev et al, 1974; Hill and Ross, 1979) 



a 



F^ h (x,Q 2 ) = 3xe 4 h -lp 



IT 



Bx(l-x)-l-x(l 



+ 



x In 



x 2 + (1 -xf+x(l -3x) 
1 + /31 



4m| 
Q 2 



x) 



— x 



4m|" 

Q 2 . 

2 8m| 



1-/3. 



(3.14) 



if the available hadronic energy squared W 2 = Q 2 (l — x)/x is larger than the production 
threshold 4m);. The NLO corrections to Eq. (3.14) are at most 20% (Laenen et al, 1994) and 
usually neglected by authors of photonic parton distributions, but a resolved contribution 
from the process 7*(Q 2 )(? — > hh 



(3.15) 



is sometimes included. Far above threshold heavy quarks are treated as light flavors with 
boundary conditions fh/j(m 2 h ) = h h {m 2 h ) = 0. 



C. Hadronic solutions 

Since the x-dependence of the boundary conditions /,/ 7 (x, Qq) cannot be calculated in 
perturbation theory, it is necessary to make theoretical assumptions. Usually one takes a 
form similar to /i/ 7 (x, Qq) = Nx a {l — x) 13 and fits the normalization N and the exponents a 
and (3 to experimental data. However, only one particular combination of photonic parton 
densities, i.e. F^(x, Q 2 ), is well constrained by experimental data from PETRA, PEP, 
TRISTAN, and LEP in the ranges 0.001 < x < 0.9 and 0.24 < Q 2 < 390 GeV 2 (Nisius, 
2000), and this combination is dominated by the up-quark density. Thus one has to rely 
further on models like Vector Meson Dominance (see the review by Bauer et al., 1978), 
which relates the photon to the p, u>, and mesons with J PC = 1~ 
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Vector Meson V 


mv/MeV 


r e y +e -/keV 


#/(4tt) 


P 


769.3 ± 0.8 


6.77 ± 0.32 


2.02 


UJ 


782.57 ± 0.12 


0.60 ± 0.02 


23.2 


<t> 


1019.417 ± 0.014 


1.297 ± 0.04 


13.96 



TABLE II. Masses, leptonic decay widths, and decay constants in the zero-width approximation 
T^ fe _ = 47ra 2 my/(3/y) for the p, u, and 4> vector mesons (Groom et al, 2000). 



|7>= E 4-\v) 



V=p,U),i 



fv ' \P P fl 



e 2 e 2 



_ g 

x (e u \uu) + e d \dd)) + — \ss) (3.16) 

U 

and allows for a successful phenomenological description of experimental data on the photo- 
production of vector mesons. A flavor SU(3)-symmetric superposition requires f p / f w = 1/3 
and f p / = —^/2/3 in fairly good agreement with experimental measurements (1/3.4 and 
— v^h3/3, see Tab. II). A coherent superposition of up- and down-quarks, which is favorable 
at large scales, is obtained with e u = 2/3 and = —1/3, while an incoherent superposition 
is obtained with e u — e^ — \. 

For safely high starting scales Ql > 2 GeV 2 , the pure VMD ansatz turns out to be 
insufficient to describe the Fl data at larger Q 2 and has to be supplemented by an additional 
hard component. For quarks it can be naturally provided by the quark-box diagram with 
four external photons, but this is unfortunately not viable for the gluon. Two solutions are 
possible. The first option is to retain a relatively large Q , Q > 1 GeV, fit the quark densities 
to F 2 7 data, and estimate the gluon input (Drees and Grassie [DG], 1985; Abramowicz, 
Charchula, and Levy [LAC], 1991; Gordon and Storrow [GS, GS96], 1992a, 1997; Hagiwara 
et al. [WHIT], 1995; Schuler and Sjostrand [SaS 2], 1995; Abramowicz, Gurvich, and Levy 
[GAL], 1998). The second option is to retain the pure VMD ansatz and start the evolution 
at a very low scale Qq ~ 0.5. ..0.7 GeV (Gluck, Reya, and Vogt [GRV], 1992b; Aurenche et 
al. [ACFGP], 1992; Aurenche, Fontannaz, and Guillet [AFG], 1994a; Schuler and Sjostrand 
[SaS 1], 1995; Gluck, Reya, and Schienbein [GRSc], 1999b). Unfortunately the parton 
distributions in the vector mesons are unknown, so that in practice one has to resort to 
those in the pseudoscalar pion. The various parameterizations are summarized in Tab. III. 

Recently, an energy-momentum sum rule for the photon has been derived (Schuler and 
Sjostrand, 1995; Frankfurt and Gurvich, 1996) 

f 1 dxx[Z(x,Q 2 ) + f gh (x,Q 2 ) + f 7h (x,Q 2 )] =1, (3.17) 

J 

which further constrains the boundary conditions. E(x, Q 2 ) = Z)i=fi 2/ g / 7 (x, Q 2 ) is the singlet 
quark distribution. Integrating the third evolution equation in Eq. (3.5) one obtains the 
solution 
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Group 


Year 


Set 


Ql 
(GeV 2 ) 


Factor. 
Scheme 


VMD 
Model 


N f 


AT. — A 

MS 

(MeV) 


DG 


1985 


1-3 


1 


LO 


- 


3-5 


400 


LAC 


1991 


1-3 


4,4,1 


LO 


- 


4 


200 


GS 


1992 


LO 


5.30 


LO 


incoherent 


5 


200 






NLO 


5.30 


MS 


incoherent 


5 


200 


GRV 


1992 


LO 


0.25 


LO 


incoherent 


5 


200 






NLO 


0.30 


DIS 7 


incoherent 


5 


200 


ACFGP 


1992 


NLO 


0.25 


MS 


coherent 


4 


200 


AFG 


1994 


NLO 


0.50 


MS 


coherent 


4 


200 


WHIT 


1995 


1-6 


4.00 


LO 


- 


4 


400 


SaS 


1995 


1D/M 


0.36 


LO 


coherent 


5 


200 






2D/M 


4.00 


LO 


coherent 


5 


200 


GS96 


1997 


LO 


3.00 


LO 


incoherent 


5 


200 






NLO 


3.00 


MS 


incoherent 


5 


200 


GAL 


1998 


LO 


4 


LO 




4 


200 


GRSc 


1999 


LO 


0.26 


LO 


coherent 


3 


204 






NLO 


0.40 


DIS 7 


coherent 


3 


299 



TABLE III. Parameterizations of the parton densities in the photon. The sets SaS 1M and 2M 
are defined by absorbing the C 7 coefficient in Eq. (3.11) into the LO quark distributions. The 
coherence of the VMD model is determined by the coefficients e u and in Eq. (3.16). Most 
parameterizations fit (except ACFGP, AFG, and GRSc) and the pion densities (except DG, 
LAC, WHIT, SaS, and GAL) and add the direct Bethe-Heitler contribution for heavy quarks in 
Eq. (3.14) (except DG, LAC, and GAL). WHIT, SaS, and GRSc also add the resolved contribution 
in Eq. (3.15). 



f lh (x,Q 2 ) = 5(1 -x) 



1 - 



a V0 J 



6tt 



(3.18) 



which can be directly inserted in Eq. (3.17) with the result 



f 

Jo 



dxx 



P Nce> m +c ) 



(3.19) 



So the quark and gluon momentum fractions rise logarithmically with Q 2 . For the sum rule 
to be of practical use, the unknown integration constant c has to be related, e.g., to the total 
hadronic cross section in e + e~ annihilation via a dispersion relation in the photon virtuality. 



Schuler and Sjostrand (1995) obtained c/tt = ^ v ^ ^/fv 
while the numerical values in Tab. II give a slightly larger value of c/n 



0.55 at Ql 



0.36 GeV 2 , 



0.61. 
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Different LO and NLO parameterizations of the parton densities in the photon are com- 
pared in Fig. 2 at a value of Q 2 = 25 GeV 2 relevant for many hard photoproduction processes. 
In LO (left) the 0(a,a s ) terms in Eq. (3.12) do not contribute and F^x^Q 2 ) is directly 
related to the quark densities. Thus the experimental constraints on F^(x, Q 2 ) lead to fairly 
good agreement of the different quark density parameterizations, particularly at large x 
(x > 0.1). Below this value, the WHIT 6 and also the similar WHIT 2-5 quark densities 
show a steep rise, which is caused by a similar rise of the gluon density. While this behavior 
is consistent with the large-a; data used in the fits, it is strongly disfavored by recent small- re 
data from LEP (Nisius, 2000). The parameterizations with a high starting scale (WHIT, 
SaS 2, and GS96) have larger quark densities at large x than those with a low starting scale 
(GRSc and SaS 1). The sets SaS 1M and SaS 2M, defined by absorbing C 7 into pointlike 
initial quark distributions even in LO, turn out similar to the sets SaS ID and SaS 2D after 
subtracting the C 7 -term. If the vector-meson states are added in a coherent fashion (GRSc 
and SaS), the d-quark input density in the photon is suppressed by a factor of four with 
respect to the -u-quark density. In this case the low-x behaviors differ, and the momentum 
sum rule Eq. (3.19) can be fulfilled. With an incoherent superposition, the -u-quark and 
rf-quark densities are identical at low x (GS96 and WHIT), and the momentum sum rule is 
easily violated. The gluon density does not enter directly in Eq. (3.12) in LO. It enters only 
through a rather weak coupling to the singlet in the evolution equations and is consequently 
badly constrained. Only a very steep rise at low x (LAC 1-2, WHIT 2-6) and a very hard 
gluon (LAC 3) can be ruled out by F^ data from LEP and jet or particle production data 
from HERA. On the other hand, the HERA data indicate that the SaS ID and the almost 
identical GRSc gluon densities may be too low (Nisius, 2000). 

The NLO parton densities are shown on the right hand side of Fig. 2 in the DIS 7 scheme. 
The AFG and GS96 quark densities have been properly transformed, while the GRSc and 
GRV densities are given directly in the DIS 7 scheme. The quark densities do therefore 
not exhibit the MS singularity at x — 1, but have shapes similar to their LO counterparts. 
Furthermore, the GRSc and AFG densities are very similar, since both groups use a coherent 
VMD ansatz at a low starting scale and relate the p-meson input to recent determinations 
of the pion densities. An incoherent superposition and older pion parameterizations have 
been used in the GRV distributions, resulting in a steeper rise for the identical -u-quark 
and ci-quark sea distributions at small x. The gluon distribution enters directly in in 
NLO (see Eq. (3.12)) and can therefore be better constrained. Again the GRV distribution 
exhibits a steeper shape than GRSc and AFG, and the GS96 result is substantially harder. 
By definition the NLO GS96 set should, like the older NLO GS set, lead to the same 
F^(x, Qq) as in LO. However, the NLO quark parameterizations turn out too small over 
the full x-range. The LO F2(x,Qq) results and the figures in the original papers (Gordon 
and Storrow, 1992a, 1997) cannot be reproduced with the available NLO parameterization, 
which is therefore not usable in its present form. 

The three usable NLO parameterizations of the photon structure function (GRSc, AFG, 
and GRV) are confronted with the world data on F^(x, Q 2 ) in Figs. 3 (small x) and 4 (large 
x). None of the parameterizations dare to describe the TPC/27 data at Q 2 = 0.24 GeV 2 , 
which may be too close to A 2 ~ 0.04 — 0.16 GeV 2 to allow for a perturbative treatment 
(see Fig. 3). The first NLO parameterization (GRV) starts at Ql = 0.30 GeV 2 . While it 
is still lower than the TPC/27 data at Q 2 = 0.38 GeV 2 , it fits the TPC/27 data rather 
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Parameterizations of the Parton Densities in the Photon 



O 

& 0.8 



o 

X 
x 



a 

x 

X 




FIG. 2. Parameterizations of the up-quark (top), down-quark (center), and gluon (bottom) den- 
sities in the photon in LO (left) and NLO (right) at Q 2 = 25 GeV 2 . The NLO AFG and GS96 
parameterizations have been transformed from the MS to the DIS 7 scheme. 
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NLO Parameterizations of the Photon Structure Function at Small x 



Q 2 = 0.24 GeV 2 

TPC/2y 
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a L3 

-▼- DELPHI 
• 



0.4 



0.2 



Q 2 = 5.0 GeV 2 



0.4 



0.2 



Q 2 = 0.38 GeV 2 
GRSc (DIS ) 
AFG (MS) 



Q 2 = 0.7, 0.71 GeV 2 




Q 2 = 10.7, 12 GeV 2 



Q 2 = 17.8 GeV 2 




Q I I I I 1 I I I I I I I 
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X 

FIG. 3. NLO parameterizations of the photon structure function at small x compared to the world 
data on (Nisius, 2000). The GRSc and GRV parameterizations are given in the DIS 7 scheme, 
while the AFG parameterization is in the MS scheme. 
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well already at Q 2 = 0.7 and 1.3 GeV 2 . The \ow-x and low-Q 2 data from OPAL and L3 
are slightly underestimated. The GRSc parameterization is evolved from Qq = 0.40 GeV 2 . 
It describes the LEP data at small x best, due to the most recent determination of the 
pion structure. AFG starts at Ql = 0.50 GeV 2 , but this parameterization should only be 
compared to data above Q 2 ~ 2 GeV 2 . The theoretical assumptions are very similar to 
GRSc, and so are the results for , AFG being generally slightly larger. At larger x and 
Q 2 (Fig. 4), the hadronic contribution decreases, while the pointlike contribution increases. 
Consequently the different parameterizations converge, and they all describe the data very 
well. The factorization scheme starts to play a role at very large x, and thus the AFG MS 
prediction differs strongly from the GRSc and GRV DIS 7 predictions. The GRSc charm 
contribution is provided only by the direct and resolved Bethe-Heitler terms in Eqs. (3.14) 
and (3.15) for Q 2 (l — x)/x > 4m 2 , while GRV and AFG have included the charm quark also 
in the massless evolution Eqs. (3.5). GRSc and AFG choose m c ~ 1.4 GeV, while GRV take 
m c — 1.5 GeV, so that the charm contribution sets on at slightly different values of x. 

Early hopes to extract the strong coupling constant from the photon structure function 
(Witten, 1977) were dashed shortly after by complications arising in higher orders (Bardeen 
and Buras, 1979; Antoniadis and Grunberg, 1983; Rossi, 1983; Gliick and Reya, 1983; Field, 
Kapusta, and Poggioli, 1987). Since then it has been widely believed that the sensitivity of 
the photon structure function to the strong coupling is low. Most of the parameterizations 
use an approximate NLO solution of the renormalization group equation 



with the one- and two-loop /5-functions (3® — 11 — 2Nf/3 and f3\ = 102 — 38V/ /3 and the 
fundamental QCD scale parameter A^ -4 = 200 MeV. GS use the same A-value for all Nf, 
while GRSc take a s (Mz) = 0.114 and solve the renormalization group equation numerically. 
However, a new analysis shows that the now final PETRA, TRISTAN, and LEP data lead 
to a competitive determination of a s from a single-parameter pointlike fit at large x and Q 2 
with a s (mz) = 0.1183 ±O.OO5O(exp.)lo;o02g(theor.) and from a five-parameter full (pointlike 
and hadronic) fit at all x and Q 2 with a s {m z ) = 0.1198±0.0028(exp.)±^4^theor.) (Albino, 
Klasen, and Soldner-Rembold, 2002). 

The hadronic structure of the photon is currently much less well known than that of the 
proton, and current parameterizations rest largely on similar theoretical assumptions like 
Vector Meson Dominance. To improve the situation, the light quark and gluon densities have 
to be disentangled from various observables in jet, hadron, and prompt photon production 
(see the following Sections), and the heavy quark contributions have to be consistently 
resummed according to variable flavor number schemes. 



a s (Q 2 ) 




(3.20) 
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NLO Parameterizations of the Photon Structure Function at Large x 
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FIG. 4. Same as Fig. 3, but now for large x. 
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IV. JET PRODUCTION 



Due to the confinement properties of QCD, only hadrons, and not partons, are observable 
as asymptotic states. The hadronization process happens at low scales, where the coupling 
a s is large, and is therefore not calculable in perturbation theory. It can, however, be 
modeled by combining several partons or hadrons moving in the same direction into jets. 
At sufficiently large transverse energies, the production cross section can then be calculated 
perturbatively and jets can be related to partons. 

Jets were first observed in event deformations at the e + e~ colliders SPEAR (Hanson et 
al, 1975) and DORIS (Berger et al, 1978), where they helped to identify the quark spin. 
At PETRA, they could then be seen even by the naked eye. They established the existence 
of the gluon (Bartel et al, 1980; Brandelik et al, 1980; Berger et al, 1980; Behrend et al, 
1982) and were also observed in photon-photon collisions (Bartel et al, 1981; Brandelik et 
al, 1981; Althoff et al, 1984; Berger et al, 1984a, 1985, 1987; Behrend et al, 1991). Jet 
production in photon-hadron collisions was subsequently observed at Fermilab (Adams et 
al, 1994) and at HERA (Abt et al, 1993; Derrick et al, 1995a). In the following, various 
jet definitions and jet cross sections will be discussed. 

A. Cone and cluster algorithms 

In Sec. Ill we showed, how photons can split into massless quark-antiquark pairs and 
exhibit a collinear singularity. Collinear divergencies arise quite generally in scattering 
processes, whenever a massless parton splits into two. In addition, infrared divergencies 
are generated by the emission of soft particles. Within the context of QED, it could be 
proven that soft (Bloch and Nordsieck, 1937) and collinear (Kinoshita, 1962) divergencies 
cancel between real and virtual emission in total cross sections. Lee and Nauenberg (1964) 
demonstrated, that the cancellation happens also for differential cross sections, if they are 
summed over degenerate initial and final states. In QED the collinear divergency is regulated 
by the electron mass, but in QCD it occurs even for non-zero quark masses, since one massless 
gluon can split into two. Nevertheless, infrared and collinear singularities are expected 
to cancel for all leading-twist QCD observables. A rigorous proof of the Kinoshita-Lee- 
Nauenberg (KLN) theorem does not exist for QCD, but in practical applications counter- 
examples have never been found. 

Shortly after the first experimental observation of jets, Sterman and Weinberg (1977) 
realized, that beyond LO there can be no one-to-one correspondence between jets and par- 
tons, but that jet definitions must be infrared safe and satisfy the KLN theorem order by 
order in perturbation theory. In practice, jet definitions should be simple to implement in 
experimental analyses and theoretical calculations and be insensitive to hadronization ef- 
fects. In the Sterman- Weinberg definition, a final state in e + e~ annihilation is classified as 
two-jet like, if all but an energy e is contained in a pair of cones of half-angle 5, while all 
other events are classified as three-jet like. This cone definition is infrared safe, but it is 
not well-suited for events with more than three jets, and it separates the phase space in an 
inefficient way. A better solution is provided by clustering algorithms such as the one used 
by the JADE collaboration (Bethke et al, 1988), where two particles % and j are combined 
into a jet and their four-momenta p it j are added, if their invariant mass 
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(Pi + Vjf = 2EiEj(l - cos %) < yS (4.1) 

is smaller than some fixed fraction y of the overall center-of-mass energy S. Theoretically the 
JADE algorithm is not well suited to resum logarithms of 0(a™ In 2 ™ y), which can become 
large for small values of y. This is easier in the kr or Durham algorithm (Catani et ai, 
1991), where two particles are combined, if their energies Eij and the angle % between 
them satisfy the condition 

2 min(£ 2 , £ 2 )(1 - cos %) < yS. (4.2) 

In photon and hadron collisions, the partonic system is boosted along the beam axis, and 
the remnants of the beams must be separated from the hard jets. It is therefore convenient 
to define outgoing particles i in the transverse plane with four-momenta 

Pi = (E i ,p x .,p y .,p g .) (4.3) 
= (m Ti cosh y h p Ti cos fa, p Tl sin fa, m Ti sinh yi) , 



transverse masses = \jp\ + m i , transverse momenta p^, azimuthal angles fa, and 
rapidities 




Ej + Pz 
Ei -p z 



(4.4) 



Experimentally the jet energies Ei and scattering angles 6i are measured, and the are 
replaced with E Ti = i^sinflj. The rapidities yi are additive under boosts along the z- 
direction. For massless particles they coincide with the pseudorapidities rji = — ln[tan(6 l i/2)]. 
According to the standardization of the 1990 Snowmass meeting (Huth et ai, 1990), particles 
i are added to a jet cone J, if they have a distance 

Ri = ^(^-Vjy + ^-^j) 2 < R (4.5) 

from the cone center. R = 0.7... 1 is the jet radius in the rj — <fi plane and 

ETj = E E Ti , VJ = p7~ E E TiVi, 
Ri<R T J Ri<R 

and 0j = -J- E tA (4-6) 

^Tj Ri<R 

define the jet transverse energy and axis. In NLO QCD two partons, which may be separated 
by as much as 2R, can be combined, if they have equal transverse energies. In this case, one 
could count the individual jets in addition to the combined jet, but this double-counting 
should be avoided (Ellis, Kunszt, and Soper, 1989b). The broad combined jet is difficult to 
find experimentally, since it does not have a seed in its center. This led Ellis, Kunszt, and 
Soper (1992) to propose an additional parameter -R sep , which restricts the distance 



Rij = y(vi - Vj) 2 + (fa - <A?) 2 



< mm 



Et ~\~ Et 



max(E T .,E T )' 



(4.7) 
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of two partons i and j. However, the value of R scp = R...2R must be fitted to different 
experimental implementations of the Snowmass algorithm (Butterworth et ai, 1996; Klasen 
and Kramer, 1997b). Furthermore, studies of three jet hadroproduction (Giele and Kilgore, 
1997) and jet shapes (Seymour, 1998) revealed, that the Snowmass cone algorithm is not 
infrared safe in higher orders, unless the centers of jet pairs are also considered as seeds 
(Akers et ai, 1994). 

The deficiencies of the cone algorithm are remedied in the longitudinally invariant kx 
clustering algorithm (Catani, Dokshitzer, and Webber, 1992; Catani et ai, 1993; Ellis and 
Soper, 1993), where one uses only the combination criterion < R sep with the .Bp- weighted 
recombination scheme as in Eq. (4.6). The results obtained with the cluster algorithm in 
hadron-hadron collisions become similar to those obtained with the cone algorithm, if one 
chooses -R S ep — 1.35.R (Ellis and Soper, 1993). This was also found to be the optimal value 
of _R sep for the cone algorithm in photoproduction (Butterworth et ai, 1996; Klasen and 
Kramer, 1997b). 



B. Single jets 

Since photons couple to charged quarks, they can participate directly in hard scattering 
processes like jet production. On the other hand, they can also exhibit a partonic struc- 
ture (see Sec. Ill) and interact indirectly through their constituents (quarks and gluons). 
In photon-photon scattering, this leads to the direct, single- resolved, and double-resolved 
processes shown in Fig. 5, where the hard central jets are produced in association with 
zero, one, or two forward-going photon remnant jets. In photon-hadron scattering, the 
single- and double-resolved photon processes described above correspond to the direct and 
resolved photon processes in Fig. 6, which involve a proton remnant jet and, for resolved 
photoproduction, also a photon remnant jet. 

The LO partonic cross section for the scattering of two massless partons a and b into 
two massless partons 1 and 2 is 



dt 2s T(l - e) V tu J 8ns S a S b C a C b 



da B 1 1 /4tts\ £ 1 gl b 



M B \ 2 (4.8) 



in n = 4 — 2e dimensions, s = (p a + p b ) 2 , t = (p a — pi) 2 , and u = (p a — P2) 2 are the Lorentz- 
invariant Mandelstam variables, T(x) is the Euler T-function, and g 2 b are the squared cou- 
pling constants (4nae 2 for photons and Ana s {ji 2 ) for quarks and gluons). S a>b are the spin 
degrees of freedom of the initial particles (2 for quarks and n — 2 for photons and gluons) 
and C a)b are the initial color degrees of freedom (1 for photons, Nc for quarks, and Nq — 1 
for gluons). While one can safely set n = 4 or e = in Eq. (4.8) and in the squared matrix 
element |.M B | 2 in LO calculations (Brodsky et ai, 1978, 1979; Owens, 1980; Fontannaz et 
ai, 1980; Baer, Ohnemus, and Owens, 1989a), the 0(e) terms are needed in NLO in connec- 
tion with ultraviolet, infrared, and collinear divergencies (see below). In jet cross sections, 
the partonic masses p 2 = m 2 <C E\ can usually be neglected, and the Mandelstam variables 

s = Ax a x b E A E B , t = -2x a E A E T e~ v , 

and u = -2x b E B E T e r] (4.9) 
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FIG. 8. Born diagrams for single-resolved jg scattering. 



become simple functions of Et, rj, Ea,b, and x a ^. 

The massless Born diagrams for the direct process 77 — > qq, the single-resolved resolved 
process jg — * gg, and the double-resolved processes qq' — > qq', qq — > qq, qq — > (7(7, and 
yfl 1 ~~ *■ ^5 are shown in Figs. 7, 8, and 9, and the corresponding squared matrix elements 
|.M B | 2 are summarized in Tab. IV. All other Born diagrams and squared matrix elements 
can be obtained by crossing particle lines and Mandelstam variables and including a factor 
of (—1) for every crossed fermion line. 

In NLO, the partonic cross section 



da 
dt 



da B da v da F da 1 
+ —r- + —r- + 



dt 



dt 



dt 



dt 



(4.10) 



receives additional contributions from virtual loop corrections (V) and real emission correc- 
tions in the final state (F) and initial state (I) with 



da' 



a 2 
9a,b 



dt 



X 



/47TS\ £ 1 

e) V tu J 8ns S a S},C a Cb 
" \M V ^\ 2 . 




(4.11) 



Virtual corrections arise in NLO by dressing all constituents of the LO diagrams in 
Figs. 7, 8, and 9 with one-loop particle exchanges. Fig. 10 shows the self-energy diagrams 
for quarks (top) and gluons (bottom), which are factorizable and thus independent of the 
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FIG. 9. Born diagrams for double-resolved qq', qq, qq, and gg scattering. 



Process 


LO matrix element squared A4 B 2 


77 QQ 
19 99 
gg' -> gg' 

gg -> gg 

QQ^ 99 
99 -»• 55 


( 

i 


lN c (l-e) [(1 - e ) (a + ±) _ 2e] 
liVcC F - e) 

'|M B |^,,(a,t,«) + IM^^s.M) 
-8C7Hl-e)(S+e)]/2! 

[4C F (1 e)( 2iV ^ ^)(t 2 + ^ 2 ^ 2 )]/2! 
32JV& C F (1 - e) 2 (3 - $ - f - §)] /2! 



TABLE IV. LO squared matrix elements |.M B | 2 for massless 2^2 parton processes involving 
two, one, and no photons in n = 4 — 2e dimensions. 



27 




FIG. 10. QCD quark (top) and gluon (bottom) self-energy diagrams. Faddeev-Popov ghosts are 
depicted as dashed lines. 



scattering process and have to be inserted in external and internal quark and gluon lines. 
The ghost self-energy is obtained by replacing quark lines with ghost lines in the upper 
diagram. While photon self-energies arise also at one loop, they are of 0(a), not of 0(a s ), 
and thus not included in NLO QCD calculations. Vertex corrections are shown in Fig. 11 
for the quark-gluon (top), triple-gluon (center), and four-gluon vertices (bottom). For the 
latter two, additional diagrams are generated by permutations of the participating particles. 
The corrections to the ghost-gluon vertex are of similar type as those for the quark-gluon 
vertex, and the photon-quark vertex correction is obtained from the second diagram in the 
first row of Fig. 11 by replacing the outgoing gluon with a photon. Vertex corrections 
depend on the four-momentum flow at the vertex and thus factorize only at the amplitude, 
and not the squared amplitude, level. A third type of loop diagrams is generated by box 
diagrams, which depend on the types and momenta of all parton legs in a 2 — > 2 Born 
process and factorize also only at the amplitude level. Typical examples are the boxes 
generated by the four-gluon vertex correction in the last row of Fig. 11. All loop diagrams 
with a closed fermion line receive a factor (—1) and those with identical particles (gluons) 
a statistical factor of 1/2. | A^ v | 2 in Eq. (4.11) denotes the interference terms of Born 
and virtual one-loop matrix elements, which have been integrated over the loop momentum 
I with a measure d n //(27r)™. In this case only the real parts of the loop integrals have 
to be kept. In loop diagrams with internal fermion lines, the integration momentum I 
appears not only in the propagator denominator, but also in the tensor-valued numerator. 
In 2 — > 2 scattering these integrals can involve between one and four vertices. Tensor 
integrals can be reduced to scalar integrals by exploiting Lorentz invariance and extending 
the well-known reduction procedure (Passarino and Veltman, 1979) to n dimensions. The 
coefficients of the scalar integrals are then finite and have to be kept up to 0{e) or 0{e 2 ). 
All divergences are contained in the scalar integrals. They are of three different types. 
Ultraviolet 1/e poles appear at the upper boundary (oo) of the energy integration in the 
one- and two-vertex functions. They are removed, together with the universal finite terms 
— 7# + ln(4-7r), by renormalizing the bare strong coupling g s , wave functions, and masses rhi 
(where applicable) in the MS scheme (Bardeen et ai, 1978). Infrared poles arise at the lower 
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FIG. 11. QCD corrections to the quark-gluon vertex (top), triple-gluon vertex (center), and four- 
gluon vertex (bottom). 




FIG. 12. Real QCD corrections through 1 — >• 2 parton splitting. 



boundary (0) of the energy integration, when massless particles are exchanged between real 
particles, and collinear poles occur in the splitting of massless particles into two massless 
collinear particles. Infrared and collinear poles and double poles appear in the two-, three-, 
and four-vertex functions and in the derivatives of two-point functions, which play a role 
in external self-energy corrections. These poles are particularly cumbersome in massless 
QCD calculations. In practice, the scalar integrals can be calculated either with Feynman 
parameters ('t Hooft and Veltman, 1979) or with Cutkosky cutting rules and dispersion 
relations ('t Hooft and Veltman, 1973). Analytical continuation of products of logarithms 
can lead to numerically large terms of 7r 2 . The complete one-loop QCD corrections have been 
calculated in dimensional regularization for massless parton scattering processes involving 
two (Aurenche et al, 1984a), one (Aurenche et al, 1987), or no photons (Ellis and Sexton, 
1986). 

Real corrections arise in NLO by the splitting of one particle into two as shown in Fig. 
12 for massless QCD partons. These diagrams have to be attached at all possible places to 
the LO diagrams in Figs. 7, 8, and 9 and lead to 2 — > 3 processes with different kinematics. 
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The first three diagrams could also be applied to photons, but they would constitute 0(a) 
corrections and should not be included at 0(a s ). Only the 7 — > qq diagram plays a special 
role, since it provides a link at NLO between direct and resolved photoproduction (see 
Sec. III). Analytical expressions for the massless 2-^3 scattering matrix elements have 
been computed in n dimensions for photon-photon (Aurenche et al, 1984a), photon-parton 
(Aurenche et al., 1987), and parton-parton (Ellis and Sexton, 1986) scattering. 

For single-particle production (see Sees. V and VI), it is possible to integrate the squared 
2^3 matrix elements |.M F,/ | 2 analytically over the solid angle 

fdn= r d0i sin 1 " 2 " 6 l r &6 2 sin~ 2e 6 2 (4.12) 
J Jo Jo 

of the two unresolved final (F) or initial (I) state particles in a frame, where they have no 
net three-momentum. However, the integration over the invariant mass s^- = (pi + Pj) 2 of 
the two unobserved partons with four-momenta pi and pj has to be performed numerically 
(Ellis et al, 1980). These integrations are implied in the cross section expression of Eq. 
(4.11). The method described above can be applied to direct (Aurenche et al, 1994a), 
single- (Aurenche, Fontannaz, and Guillet, 1994b) and double-resolved jet photoproduction 
(Aversa et al, 1989, 1990), if the phase space in Eq. (4.12) is restricted to allow for an 
observed jet with a small cone of angular size 8 similar to the Sterman- Weinberg criterion. 
For comparisons with experimental data, finite cone size corrections have to be implemented 
by integrating the 2 — » 3 matrix elements numerically in the finite region between the small 
technical cone of size S and the larger experimental cone of size R (Gordon and Storrow, 
1992b; Greco and Vicini, 1994; Aversa et al, 1991; Kramer and Salesch, 1994a; Aurenche et 
al., 1994a; Bodeker, Kramer, and Salesch, 1994; Aurenche, Fontannaz, and Guillet, 1994b; 
Klasen, Kramer, and Salesch, 1995). The double- resolved contributions in these calculations 
were all based on the same results by Aversa et al. (1989, 1990). With the exception of 
the first two calculations, they were found to agree numerically with later calculations using 
different methods (Klasen, 1996b; Aurenche et al., 2000). Similar methods have been applied 
quite early in single-resolved calculations (Baer, Ohnemus, and Owens, 1989b; Gordon and 
Storrow, 1992b; Bodeker, 1992a, 1992b, 1993) and will be described in detail in Sec. IV.C. 

If the radiated gluons in Fig. 12 become soft and/or two of the three particles become 
collinear, singularities appear in the propagators of the intermediate particles and have to be 
extracted analytically as \ je or 1/e 2 (double-)poles. The interference of different amplitudes 
can furthermore involve multiple singularities that have to be separated by partial fraction- 
ing. According to the KLN theorem, all soft and final state collinear singularities from the 
real corrections have to cancel against corresponding virtual singularities in inclusive jet 
cross sections. Only universal collinear singularities 

l-M 7 |^i 23 = -- £ PcUx)\M B \ 2 cb ^ 12 + O(e ) (4.13) 

remain, which correspond to the parton splittings in Fig. 12 attached to the initial state 
of the Born diagrams. They are proportional to the Altarelli-Parisi splitting functions in 
Eq. (3.6) and are absorbed, together with scheme- dependent finite terms, into the parton 
densities in the proton or in the photon (see Sec. III). For the production of jets in direct, 
single-resolved, and double-resolved photon-photon collisions, the soft and collinear singular 
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Process 


Color 
Factor 


Correc- 
tion 


Singular Parts of 

\ M V,F,I\2/\ M B\2 


77 QQ 


c F 


Virtual 
Final 




-£-I(3-2Z(f)) 
+£ + I(3-2Z(t))" 



TABLE V. Cancellation of soft and collinear singularities from virtual and final state NLO cor- 
rections for direct photon scattering. 



parts of the NLO squared matrix elements |.A/f ^'^l 2 are collected in Tabs. V, VI, and VII. In 
all three cases the sums are finite, since the factorizable initial state singularities have already 
been subtracted. The arguments of the logarithms l(x) = \n\x/Q 2 \ have been normalized 
to the arbitrary scale Q 2 , which appears also in Eq. (4.11) and drops out in the total result. 



Process 


Color 
Factor 


Correc- 
tion 


Singular Parts of 
\M V ^\ 2 /\M B \ 2 


1Q 9Q 


c F 


Virtual 

Final 

Initial 




-£-I(3-2Z(f))j 
+^ + ^(3-2/(t))" 
"+£ + £(3-2^))" 




N c 


Virtual 

Final 

Initial 




~ h ( t ~ 2l ( s 

+ h H(*A) ~ 1(U) 


It) ~ 2l{u))\ 

t)-l(u))' 

} 


N f 


Virtual 
Final 


1 

3e 



TABLE VI. Cancellation of soft and collinear singularities from virtual, final state, and initial 
state NLO corrections for single-resolved photon scattering and different color factors. 



The photoproduction cross section for inclusive jets 

TrTi = S / dx a x a f a/A (x a , M 2 )x b f b/B (x b , Ml) 
AEaEt 



x 



da 



2x a E A - E T e^ dt 



(4.14) 



depends on the transverse energy E T and pseudorapidity r\ of the observed jet. In Eq. (4.14), 
f] is defined in the laboratory frame with particle A travelling in the positive z-direction. It is 
related to the pseudorapidity in the hadronic center-of-mass frame rj* = rj + 1/2 \yl(Eb/ E A ) 
by a simple boost. From the final state variables E? and rj and from the energies E A 
and E B of the initial leptons or hadrons, the four-momenta of the incoming partons p a ^ b = 
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Process 


Color 
Factor 


Correc- 
tion 


Singular Parts 
|^W| 2 /|M B 


of 

2 




qq' -> qq 1 


C F 


Virtual 

Final 

Initial 




-^-i(6 + 8/(s/«)-4/(t)) 
"+Jr + ±(3 + 4/(s/u)-2/(i))~ 
+£ + §(3 + 4l(s/u)-2i(t))" 






Virtual 

Final 

Initial 




+±(4Z(s) - 2Z(«) - 2l(t) 
'-1(21(8) - l(u) - l(t)) 
-1(21(8) - l(u) - l(t)) 


)J 


qq -> 


Cf 


Virtual 

Final 

Initial 




-£-i(6 + 4J(s/t)-4Z(u)) 
+|r + i(3 + 2i(s/t)-2Z(u))' 
"+£ + ±(3 + 2i(s/t)-2Z(u))' 




7V C 


Virtual 

Final 

Initial 




+l(2l(s)-l(t)-l(u)) 
' -1(21(8) -l(t)-l(u)) 
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TABLE VII. Cancellation of soft and collinear singularities from virtual, final state, and initial 



EA,B(x a ,b, 0, 0, ±x 0j &) cannot be fully reconstructed. Instead the partonic cross section da/dt 
is integrated over x a with 



2E A E B - E A E T e~^ ' 



x = x a E AE T e v 

Xb 2x a E A E B - E B E T e^ 1 ' ' 

and convolved with the parton densities f a / A (x a , M 2 ) and fb/B(xb, M 2 ) at the factorization 
scales M a ^. The density of partons i in initial leptons / is given by the convolution 

f ifl (x,M 2 ) = \f lh {M 2 ) ® / 7/J ](a;) (4.16) 

of the photon spectrum f 7 /i(x) as described in Sec. II with the parton densities in the photon 
fi/ y (x, M 2 ) (see Sec. Ill), so that the two can also not be uniquely disentangled. 

Due to the truncation of the perturbative series, fixed order predictions depend on the 
scales, at which the strong coupling a s (fi) is renormalized and at which the collinear initial 
state singularities are factorized into the parton densities fi/ a (x, M 2 ). In LO, the logarith- 
mic scale dependence is generally quite strong. It is reduced in NLO through the explicit 
logarithmic dependences in the virtual and real initial state corrections 

M V \ 2 ab , 12 lnf^l) M' ! :„, 1 {] (1.17) 




lA^iaa = 1^ gjfj \M B \ 2 cb ^ 12 P c Ux)+ - • (4-18) 

For the renormalization scale //, this is demonstrated in Fig. 13, where the total single- 
jet photoproduction cross section Eq. (4.14) has been calculated using HERA kinematics 
(E p = 820 GeV and E e = 26.7 GeV). Unless stated otherwise, NLO parton densities in the 
MS scheme for the proton (Lai et al, 1995) and photon (Gliick, Reya, and Vogt, 1992b) 

Nt =4 

have been used with A^ = 239 MeV. For a particular choice of //, the NLO corrections 
vanish and the perturbative series converges apparently very fast (Grunberg, 1980). If a 
one-loop a s formula is used in the LO calculation, this happens at /i ~ M p = M 1 = Et, 
while for two-loop a s it happens at fi ~ Et/4. At this scale, the NLO curve exhibits also 
a local maximum and has minimal sensitivity to scale variations (Stevenson, 1981a, 1981b). 
However, the dependences on M p and M 7 also have to be considered and can spoil this 
simple picture. If all scales are varied simultaneously, the Principles of Fastest Convergence 
and Minimal Sensitivity lead indeed to unreasonably low scale choices (Klasen, 1996b). The 
situation becomes even more complicated, if more than one scale is involved, as is the case 
in heavy quark or DIS jet production. 

In photoproduction, the dependence on the photon factorization scale M 7 is of particu- 
lar interest, since this scale defines the separation of direct and resolved processes. While 
the direct LO cross section is manifestly independent of M 7 , the resolved LO cross section 
depends logarithmically on M 7 through the LO parton densities in the photon. This depen- 
dence is canceled almost exactly by the explicit logarithmic dependence of the NLO direct 
term, as can be seen in Fig. 14 (Bodeker, Kramer, and Salesch, 1994; Klasen, 1996b; Klasen, 
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FIG. 13. Renormalization scale dependence of the single-jet cross section dV/dEydr/ at HERA 
with Et = 20 GeV and rj = 1 in LO (dotted and dashed curves) and NLO (full curve). 



Kleinwort, and Kramer, 1998). However, the resolved cross section rises faster in NLO than 
in LO (Klasen, 1996b). This reintroduces a (weaker) dependence in the total NLO result, 
which would eventually be canceled by the direct NNLO contribution. 

The HERA experiments HI and ZEUS select photoproduction events by tagging and anti- 
tagging the scattered electron, respectively, which limits the maximal squared virtuality of 
the exchanged photon to Q 2 < 0.01.. .1 and 1...4 GeV 2 . The energy fraction y transferred to 
the photon is determined from the scattered electron energy E' e with y — 1 — E' e /E e or with 
the Jacquet-Blondel method 

2/=^E(^-^) ( 4 - 19 ) 

from the energies Ei and longitudinal momenta p Zi of the final hadrons (Amaldi et al, 1979). 
Single-jet production has been analyzed with the cone algorithm and R = 1 by HI (Abt 
et al, 1993; Aid et al, 1996a) and ZEUS (Derrick et al, 1995a; Breitweg et al, 1998a) 
and has been compared to various NLO calculations (Aurenche, Fontannaz, Guillet, 1994b; 
Klasen, Kramer, and Salesch, 1995; Klasen, 1996a, 1996b, 1997a; Klasen, Kleinwort, and 
Kramer, 1998; Harris and Owens, 1997, 1998). In these early comparisons at rather low 
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FIG. 14. Photon factorization scale dependence of the single-jet cross section d 2 cr/d-E/rd?7 at HERA 
with Et = 20 GeV and rf = 1. The LO direct curve (dotted) is independent of M 7 as is the sum 
(full curve) of NLO direct (dashed) and LO resolved (dot-dashed) contributions. 



transverse energies of ~ 10 GeV, the theory overestimated the data in the backward region, 
since no hadronization corrections had been applied, while they underestimated the data 
in the forward region, triggering speculations about additional soft interactions between 
the hard jets and the proton remnant. In more recent analyses (Breitweg et al, 1998b; 
Adloff et al, 2001a), the invariant kx algorithm has been used. The measured HI cross 
section in Fig. 15 now extends out to < 75 GeV for large y. It falls by three orders of 
magnitude and agrees remarkably well with the NLO prediction (Frixione and Ridolfi, 1997). 
The hadronization corrections $hadr. are m uch smaller than the NLO corrections, and the 
remaining theoretical uncertainty from varying the scale /j, = M p = M 7 = E T /2 by a factor of 
two up and down is modest. Different photon densities give slightly different normalizations, 
but the experimental uncertainty is still too large to draw definite conclusions. Single-jet 
pseudorapidity distributions are shown in Fig. 16. The maximum of the cross section is 
shifted in the electron direction (low rj) when increasing y and lowering E T . The HI data are 
in general well described by the NLO prediction, perhaps best with GRV photon densities, 
but the experimental and theoretical (not shown) uncertainties are still quite large. Direct 
photoproduction dominates at large Et and small rj, while resolved processes are more 
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FIG. 15. Transverse energy dependence of the total single-jet photoproduction cross section com- 
pared to preliminary HI data (Adloff et al, 2001a). 
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FIG. 16. Pseudorapidity dependence of the total single-jet photoproduction cross section compared 
to preliminary HI data (Adloff et al, 2001a). 
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important at smaller Et and larger rj. 

A good measure for the width of a produced jet is the jet shape 



<R) = 1 _ I<iE T ,l„E T ,l^o/(,lE T d,) 

which is the average fraction of a jet's transverse energy that lies inside a concentric inner 
cone with radius r < R. It has been calculated in resolved (Kramer and Salesch, 1993, 1994b; 
Klasen and Kramer, 1997b; Klasen, 1997a) and direct (Klasen and Kramer, 1997b; Klasen, 
1997a) photoproduction by integrating the ^-weighted NLO 2^3 cross section in the 
non-singular region between r and R. In Fig. 17 jet shapes for single-jet photoproduction 
are compared to ZEUS data (Breitweg et al, 1998c) for different transverse energy and 
pseudorapidity intervals. It is obvious that the data and the NLO predictions rise more 
steeply for large values of Et and small values of i], where the jets become narrower. Using 
LO Monte Carlo predictions, the ZEUS collaboration has recently been able to associate 
thick and thin jets and with gluon and quark jets (Breitweg et al, 2000a). 

Inclusive jet production in photon-photon scattering has been measured at the e + e~ 
collider TRISTAN (y/S = 58 GeV) by the TOPAZ (Hayashii et al, 1993) and AMY col- 
laborations (Kim et al, 1994), where scattered electrons were anti-tagged at small angles 
(#max = 227 and 56 mrad). OPAL (Ackerstaff et al, 1997) performed a measurement at 
LEP1.5 (\/S = 133 GeV) with a maximal electron scattering angle of 33 mrad. All exper- 
iments used the cone algorithm and R — 1. As shown in Fig. 18, the AMY data are well 
described by the sum of NLO direct, single-, and double-resolved processes (Kleinwort and 
Kramer, 1996a) with GRV photon densities. The same observation is made in NLO compar- 
isons with the TOPAZ (Aurenche et al, 1994a; Kleinwort and Kramer, 1996a) and OPAL 
data (Klasen, Kleinwort, and Kramer, 1998). At comparable values of xt = 2Et/\/S, the jet 
width measured in photon-photon collisions is approximately the same as in photon-hadron 
collisions (Adachi et al, 1999). 



C. Dijets 

NLO calculations for dijet photoproduction differ significantly from those for single jets 
with respect to the real corrections. While the phase space of the two unobserved partons 
could be integrated analytically for single jets (see Eq. (4.12)), it is now restricted by the 
definition of the second observed jet. The third unobserved parton momentum p 3 still has 
to be integrated out. 

In the phase-space slicing method, this is done analytically in the soft and collinear 
regions, which can be defined by limiting the invariant masses of two unresolved partons 
Si3 = (Pi + P3) 2 < ys (Gutbrod, Kramer, and Schierholz, 1984; Kramer and Lampe, 1989) 
or by applying separate soft and collinear cut-offs on E 3 < ey/s/2 and 6 i3 < 5 (Fabricius et 
al, 1981; Gutbrod, Kramer, and Schierholz, 1984), on E 3 < 5 s ^fs/2 and s i3 < 5 c s (Baer, 
Ohnemus, and Owens, 1989b; Harris and Owens, 2001), or on Et 3 < PTm and Ri 3 < R c 
(Aurenche et al, 2000). Outside the singular regions, the real corrections are integrated 
numerically, so that the cross sections become independent of the technical cut-offs and an 
experimental jet definition can be implemented. While the y-cut method has been applied 
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FIG. 17. Jet shapes p(r) for single-jet photo-production integrated over — 1 < r/ < 2 and four 
different regions of Et (a)-(d) and integrated over Et > 14 GeV and four different regions of r/ 
(e)-(h). 1994 data from ZEUS (Breitweg et al, 1998c) are compared to NLO results using the cone 
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FIG. 18. Transverse momentum distribution of single jets produced in photon-photon collisions 
compared to AMY data (Kleinwort and Kramer, 1996a). 



to direct (Kleinwort and Kramer, 1996a, 1996b, 1997; Klasen, Kleinwort, and Kramer 1998), 
single-resolved (Klasen and Kramer, 1996a, 1996b, 1997a), and double- resolved (Klasen and 
Kramer, 1997a) dijet photoproduction, the 5 S)C (Baer, Ohnemus, and Owens, 1989b; Harris 
and Owens, 1997, 1998) andpi-m, R c methods (Aurenche et al, 2000) have only been applied 
to the single- and double-resolved cases. All these results were found to be in good agreement 
with each other (Harris, Klasen, and Vossebeld, 1999; Aurenche et al. 2000). 

In the subtraction method, p% is integrated numerically also in the soft and collinear 
regions. The 2^3 matrix elements then have to be regulated by subtracting them in 
a finite phase-space volume setting p 3 = everywhere except in the singular propagator 
denominators. The subtracted terms are then integrated analytically in the same phase- 
space volume and added to the virtual corrections to allow for an analytical cancellation 
of the soft and collinear poles (Ellis, Ross, and Terrano, 1981). Several solutions have 
been proposed to obtain numerical stability with this method (Ellis, Kunszt, and Soper, 
1989a, 1989b, 1990; Kunszt and Soper, 1992). QCD dipole factorization formulae allow 
for a straightforward computation of the soft and collinear subtraction terms (Catani and 
Seymour, 1996, 1997). The subtraction method has been applied to single- resolved (Bodeker, 
1992a, 1992b, 1993; Frixione and Ridolfi, 1997) and double-resolved (Frixione and Ridolfi, 
1997) dijet photoproduction. 
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FIG. 19. Dependences of the (a) dijet and (b) three-jet cross sections on the fastest jet scattering 
angle, normalized at cos# = 0. The ZEUS dijet (Derrick et al, 1997a) and three-jet (Breitweg 
et al, 1998d) data are well described by the QCD predictions, but not by the pure phase-space 
distributions, and they favor the single 1/t pole for massless fermion exchange over the double 1/t 2 
pole for massless boson exchange. 
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In addition to the transverse energy Et and pseudorapidity rji of the first jet, the inclusive 
dijet cross section 



^^ 2 = ^s^M)^h IB ( x ,M)% (4.2D 

depends on the pseudorapidity of the second jet 772- In LO only two jets with equal transverse 
energies can be produced, and the observed momentum fractions of the partons in the initial 
electrons or hadrons 



obs 
x a,b 



^E Ti e^/(2E AB ) (4.22) 



i=i 



equal the true momentum fractions x a ^. If the energy transfer y = E 1 /E e is known (see Eq. 
(4.19)), momentum fractions for the partons in photons x° bs = x^f/y can be deduced. In 
NLO, where a third jet can be present, the observed momentum fractions are defined by the 
sums over the two jets with highest E T , and they match the true momentum fractions only 
approximately. Furthermore, the transverse energies of the two hardest jets need no longer 
be equal to each other. Even worse, for equal Et cuts and maximal azimuthal distance 
A0 = 0i — 02 = Ti" the NLO prediction becomes sensitive to the method chosen for the 
integration of soft and collinear singularities, strongly scale dependent, and thus unreliable 
(Klasen and Kramer, 1996a; Harris and Owens, 1997; Frixione and Ridolfi, 1997). This 
sensitivity also propagates into the region of large observed momentum fractions (Aurenche 
et al, 2000). Two proposed solutions are to allow for a small difference in the theoretical 
Et cuts or for a small technical cut-off dependence, but it is preferable to cut instead on the 
average E T = (E Tl + E T2 )/2. Similarly one can define AE T = E Tl — E T2 , fj — (rji + rj 2 )/2, 
and A77 = r/x — r] 2 , which is related to the cosine of the center-of-mass scattering angle 
cos#* = tanh(A?7/2). The dijet cross section then takes the form 

d3(J 2 d3flr , An „S 

(4.23) 



dx a dx b d cos 9* s dE T dr]id7] 2 

which is particularly useful to determine the parton densities and scattering processes. 
The cos #*-dependence at HERA is shown in Fig. 19 (a) for dijet masses Myi 



(pi + P2) 2 > 47 GeV and the A>r algorithm. The ZEUS data (Derrick et al, 1997a) 
show good agreement with the QCD prediction (full curve), but clearly disagree with the 
pure phase-space distribution (dashed curve) and the Rutherford scattering form at small 
angle (1 — | cos#*|)~ 2 , which is characteristic for massless vector boson exchange in the re- 
solved processes. On the other hand, the data agree very well with the less singular form 
(1 — I cos#*|) -1 , which is typical for massless fermion exchange in the direct processes, indi- 
cating that the direct processes dominate over resolved processes in this kinematic region. 



A similar behavior is observed in three-jet cross sections with M123 = y (pi + P2 + P3) 2 > 50 
GeV for the scattering angle of the leading jet (see Fig. 19 (b)) (Breitweg et al, 1998d; 
Klasen, 1999a). NLO QCD also describes earlier ZEUS data using the cone algorithm at 
lower values of M i2 (M 12 > 23 GeV) (Derrick et al, 1996), which have been divided into 
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FIG. 20. Dependence of the dijet photoproduction cross section on the observed parton momentum 
fraction in the proton compared to preliminary HI data (Adloff et al, 2001b). For "resolved" 
photons (left), the distribution is quark- like, while for "direct" photons (right) it is gluon-like, as 
is to be expected from the contributing parton scattering diagrams. 



a natter "direct" (x° bs > 0.75) and a steeper "resolved" (x° bs < 0.75) sample (Harris and 
Owens, 1997). 

Similar analyses have recently been carried out with the algorithm by HI (Adloff 
et al, 2001b) and ZEUS (Breitweg et al, 2000b), who also measured distributions of the 
observed parton momentum fractions in the proton (only HI) and photon with asymmetric 
cuts on E Tl > 25 GeV and E T2 > 15 GeV (14 and 11 GeV in the ZEUS case). The HI x° bs 
distribution (Fig. 20) is well described by the CTEQ5M proton densities (Lai et al, 2000), 
but it is not sensitive to the poorly constrained gluon density at very small or large x p . The 
HI x° bs distribution (Fig. 21) shows a slight preference for GS96 at lower and GRV at 
larger Et (see also Adloff et al, 1998), but the experimental and theoretical uncertainties 
are still larger than the photon density variations. HI (Ahmed et al, 1995; Adloff et al., 
2000) and ZEUS (Derrick et al., 1995b; Breitweg et al., 1998e) have earlier been able to rule 
out the LO LAC1 and LAC3 photon densities in x° bs and f] distributions, which are both 
sensitive to the photon densities (Forshaw and Roberts, 1993). The ZEUS distributions in 
r\2 (Breitweg et al., 1999a) do not allow any firm conclusions, since sizeable hadronization 
corrections have not been included in the NLO predictions (Harris, Klasen, and Vossebeld, 
1999). Dijet mass distributions have been measured up to 140 (Breitweg et al., 1999a) and 
180 GeV (Adloff et al., 2001b), respectively. They show no deviations from NLO QCD. 

In photon-photon collisions, dijet transverse energy distributions have been measured 
with the cone algorithm by AMY (Kim et al, 1994), TOPAZ (Hayashii et al, 1993), and 
OPAL (Ackerstaff et al, 1997; Abbiendi et al, 1999). They are well described by NLO 
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FIG. 21. Dependence of the dijet photoproduction cross section on the observed parton momentum 
fraction in the photon compared to preliminary HI data (Adloff et al, 2001b). 



44 



QCD (Kleinwort and Kramer, 1996a, 1997; Klasen, Kleinwort, and Kramer, 1998), as are 
the OPAL rapidity distributions (Abbiendi et al, 1999). The algorithm has been used by 
ALEPH (Barate et al, 2000) and OPAL (Abbiendi et al, 2001c; Wengler, 2002) to measure 
the transverse energy and x° bs (only OPAL) distributions. The former is well described 
by NLO QCD. However, the infrared sensitivity of x° bs is reflected in Fig. 22 in the large 
fluctuation of the NLO prediction in the two highest bins. When both photons interact with 
low x° bs , multiple parton interactions (MIA) can lead to a larger observed cross section than 
predicted by NLO QCD with GRV photon densities. 

Higher center-of-mass energies may be reached at linear e + e~ colliders like TESLA, ILC 
(y/S = 500.. .1000 GeV), or CLIC (v 7 ^ = 1...3 TeV) or at a future ep collider like THERA 
(TESLA x HERA, y/S ~ 1 TeV). At linear colliders, the photoproduction cross section is 
enhanced due to beamstrahlung (see Sec. II. B) or laser backscattering (see Sec. ILC). Fig. 
23 demonstrates that tests of QCD and determinations of the photon structure could be 
largely extended: Depending on the TESLA electron beam energy (250-500 GeV) and the 
collider mode (ep or 7p), the THERA range in the average transverse energy of the two jets 
would be increased by a factor of 2-3 and the reach in x° bs by at least one order of magnitude. 
The reach in x° bs is extended by about the same amount. It would thus become possible 
to check the determinations of the gluon density in the proton obtained in deep-inelastic 
scattering experiments and to measure the gluon density in the photon down to low values 
of x (Klasen, 2001c; Wing, 2001). Similar studies have been performed for high-energy 
muon-proton collisions, although here bremsstrahlung is reduced and laser backscattering 
seems impossible (Klasen, 1997b). 



D. Three jets 

For N massless jets, one can choose 3N — 4 parameters that should span the multijet 
parameter space. They should also facilitate a simple interpretation within QCD and allow 
for a comparison of the N — 1-jet to the AT- jet cross section. In the case of = 3, the 
conventional choices are the three-jet mass M 12 3 = \J (p± + P2 + P3) 2 and four dimensionless 
parameters. The Dalitz energy fractions 

*i = TT- ( 4 - 24 ) 

specify, how the available energy is shared between the three jets. They are ordered such 
that x\ > x 2 > X3. Since x\ + £2 + ^3 = 2, only x\ and x 2 are linearly independent. The 
third and fourth parameters are the cosine of the scattering angle between the leading jet 
and the average beam direction p^v — Pa — Pb, where the incoming parton a is the one with 
the highest energy in the laboratory frame, 

Q PAVPl ( . c\r"\ 

cost*! = — — ——, (4-25) 



WavWp 



and the angle between the three-jet plane and the plane containing the leading jet and the 
beam direction, 
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FIG. 23. Dijet photoproduction at a future THERA collider as a function of the average transverse 
energy of the two jets (a) and of the observed parton momentum fraction in the photon (b). The 
thin lines in (a) show the separate contributions from the resolved (direct) processes for HERA 
and THERA with E e =250 GeV, which dominate at small (large) Ej~. 
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COS^i 



(jh x Pay)(P2 x p 3 ) 
Wi x Pav\\P2 x p 3 \ 



(4.26) 



In the soft limit, where E 3 — > and xi )2 — > 1, cos#i approaches the 2^2 center-of- 
mass scattering angle cos#*, thus relating three-jet to dijet cross sections. The similarity of 
the distributions in the dijet and three-jet scattering angle distributions has already been 
discussed (see Fig. 19). Of course, the third jet must not be too soft (or the hard jets not 
too hard) to avoid soft singularities that would have to be absorbed into the next-to-leading 
order dijet cross section. This can be achieved by a cut on x±, e.g. X\ < 0.95. However, 
since the energy of a jet is always larger than its transverse energy, a cut like £^,3 > 5 
GeV already insures the absence of soft singularities. The three jets also have to be well 
separated in phase space from each other and from the incident beams to avoid initial and 
final state collinear singularities. This is insured by a cut like | cos^il < 0.8, by cuts on the 
pseudorapidities, and by the kr jet algorithm. 
The three-jet photoproduction cross section 



has been calculated from the direct and resolved tree-level 2^3 processes by several groups 
with good mutual agreement (Baer, Ohnemus, and Owens, 1989a; Klasen, 1999a, 1999b; 
Harris, Klasen, and Vossebeld, 1999). Figure 24 shows the three-jet cross section at HERA 
as a function of the energy fractions of the leading (a) and next-to-leading (b) jet x\ and X2, 
normalized to the total cross section. The prediction from the 0(aa 2 s ) QCD matrix elements 
are rather similar to the pure phase-space distributions with constant matrix elements, but 
the ZEUS data (Breitweg et al, 1998d) slightly prefer the QCD predictions. The dependence 
on the angle between the three-jet plane and the plane containing the leading jet and the 
average beam direction is shown in Fig. 24 (c). The full QCD curve again agrees well with 
the data as do the contributions from direct photons and quarks in the photon, whereas the 
pure phase space has a completely different shape. 

Unfortunately, a full NLO calculation for three-jet photoproduction has not yet been 
performed. Such a calculation would allow, e.g., for a determination of a s from the ratio 
of three-jet to dijet cross sections and for a detailed study of the jet substructure. It could 
furthermore form the basis of a NNLO dijet calculation with reduced renormalization and 
factorization scale dependences. Higher order corrections can also be taken into account 
approximately with resummation and parton shower techniques, leading eventually to more 
precise Monte Carlo generators for use by the experimental groups. 

E. Dijets with a leading neutron 

Recently HI (Adloff et al, 2001c) and ZEUS (Breitweg et al, 2000c) have presented 
data on dijet photoproduction with a leading neutron, which is dominated by slightly off- 
shell pion exchange and can be used to constrain the parton densities in the pion. These 
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FIG. 24. Dependence of the three-jet cross section on the energy fractions x± (a) and X2 (b) and 
on the angle between the three-jet plane and the plane containing the leading jet and the average 
beam direction (c). The ZEUS data (Breitweg et al., 1998d) rule out the pure phase-space and 
pure gluon-initiated distributions. 
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are not well known, particularly at low x n and in the sea-quark and gluon sectors. The 
pion structure carries important implications for the QCD confinement mechanism and the 
realization of symmetries like isospin in nature. It is also of practical importance for the 
hadronic input to the photon structure at low scales (see Sec. III.C). 

Assuming factorization, the photoproduction cross section for two jets with a leading 
neutron 

d^(J dcr 

= E Xafa/e(x a , M 2 )x b f b/p (x b , M%)-£- (4.28) 



depends on the partonic cross section da/dt (Eq. (4.10)) and on the parton densities in the 
electron f a / e (Eq. (4.16)) and proton 

f b/p (x b ,M 2 ) = f 1 d ^~ Xn) f b/w (x n ,M 2 n )f n/p (l-x n ,t'). (4.29) 

The latter is a convolution of the parton densities in the pion f b / n with the pion flux in the 
p — > nn transition 

x [F(x n ,t')} 2 , (4.30) 

where \j{mL — t') 2 is the squared pion propagator, m n is the pion mass, and (1 — 
x n ) 1 ^ 2a ' 7t ^ t '~ m ^ accounts for the virtuality and possible reggeization of the pion. In the ZEUS 
analysis, the momentum fraction of the leading neutron x n and the momentum transfer t' 
are restricted to x n > 400 GeV/820 GeV and t' = f(pr) = —pT/x n — (l—x n )(m 2 l — x n m 2 )/x n 
with pt < x n ■ 0.66 GeV. The interaction term in the pion-nucleon Lagrangian leads to the 
numerator —t', and off- mass shell effects of higher Fock states are modeled by a form factor 

. _ J exp[b(t' - ml)}, [Exponential] 
1 > ~ \ exp[R 2 (t' - ml) /(I - x n )\. [Light Cone] 1 6 > 

Since the momentum transfer t' is small, reggeization can be neglected (a' = 0), and a 
light cone form factor with R = 0.5 GeV~ x can be chosen in good agreement with recent 
determinations (Holtmann, Speth, and Szczurek, 1996; D'Alesio and Pirner, 2000). The 
pion-nucleon coupling constant g 2 /(4n) = 14.17 can be taken from a recent extraction from 
the Goldberger-Miyazawa-Oehme sum rule (Ericson, Loiseau, and Thomas, 2000). 

In Fig. 25, LO and NLO QCD predictions (Klasen and Kramer, 2001; Klasen, 2001a, 
2001b) with SMRS3 (Martin et al, 1992) and GRS (Gliick, Reya, and Schienbein, 1999a) 
pion densities and GS96 photon densities are compared to the ZEUS measurements (Breitweg 
et al., 2000c). In the transverse energy distribution (a), the scale dependence is reduced 
considerably from LO to NLO. The normalization is sensitive to the chosen pion flux factor. 
If one includes the Regge trajectory (a' = 1 GeV -2 ) and omits the form factor, the cross 
section is reduced by 15%. The distribution in the observed momentum fraction of the 
partons in the pion (b) 
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FIG. 25. Dependence of the dijet photoproduction cross section with a leading neutron on the 
transverse jet energy (a) and on the observed parton momentum fraction in the pion (b) compared 
to ZEUS data (Breitweg et al, 2000c). The error bands show the theoretical scale uncertainty 
(medium: LO, dark: NLO) and the experimental energy scale uncertainty (light). 
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suffers from the fact, that ZEUS have applied equal cuts on E Tl ,E T2 > 6 GeV and 771,2 £ 
[—2; 2], which had to be relaxed to Et,2 > 6 GeV— A in the NLO calculation. The best fit is 
obtained with GS96 photon densities and SMRS3 pion densities, which have the lowest quark 
and gluon distributions, respectively, but hadronization corrections may play an important 
role at small x n . 
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V. HADRON PRODUCTION 



In the last Section it was demonstrated, that jets provide an intuitive link with the 
partonic scattering processes, yield large rates in experimental measurements, and are well 
suited to determine the structure of the incoming photons, protons, or pions. However, they 
do not provide detailed information about the final state hadronization, which can be better 
studied in the production of light and heavy hadrons and quark-antiquark bound states 
(quarkonia). 

Early measurements of inclusive light hadron production were performed in fixed target 
collisions by the CERN NA14 experiment (Auge et ai, 1986; Barate et ai, 1986) and 
in photon-photon collisions at DESY PETRA (Berger et ai, 1979; Brandelik et ai, 1981). 
Heavy mesons were studied by the CERN WA4 (Roudeau et ai, 1980) and NA14/2 (Alvarez 
et al, 1992, 1993), Fermilab E687 (Frabetti et al, 1993, 1996; Moroni et al, 1994) and E691 
(Anjos et al, 1989, 1990), DESY JADE (Bartel et ai, 1987) and KEK AMY (Aso et ai, 
1995; Takashimizu et al, 1996), TOPAZ (Enomoto et ai, 1994a, 1994b), and VENUS 
experiments (Uehara et ai, 1994, Ohyama et ai, 1997). 

In this Section, the fragmentation of partons to light and heavy hadrons and quarkonia 
and their production in photon-hadron and photon-photon collisions will be discussed. 

A. Fragmentation 

As mentioned in Sec. IV.A, the KLN theorem guarantees, that the soft and collinear sin- 
gularities generated by real particle emission beyond the LO cancel against those generated 
by virtual particle exchanges after summation over degenerate final states. In semi-inclusive 
final states, which contain identified particles, this cancellation is incomplete. Collinear 
1/e-poles multiplying the time-like partonic splitting functions Pj<-i(x) in the transition 
functions 



remain. As in the case of parton distributions (see Sec. Ill), the factorization theorem (Ellis 
et ai, 1979; Baier and Fey, 1979; Altarelli et ai, 1979; Furmanski and Petronzio, 1980, 
1982) allows to absorb these singularities in the MS-scheme into renormalized fragmentation 
functions 



where x is the longitudinal momentum fraction of the hadron H in the parton i and 
DH/i(x), Dn/j(x) are the bare fragmentation functions. In this way one obtains important, 




(5.1) 



D H/i (x, Mj) = D H/i (x) + T^i(M 2 f ) <g> D H/j (x), 



(5.2) 
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yet incomplete, perturbative information about the confinement of unobservable quarks and 
gluons into observable hadrons. The fragmentation functions satisfy the sum rules 



£ f dxxD H/i (x,M 2 ) 



1, 



(5.3) 



i.e. the momentum of the parton % must be conserved after hadronization into all avail- 
able hadrons H. Fragmentation functions constitute only the leading-twist contributions 
to semi-inclusive hadron production in the operator product expansion. Higher- twist, non- 
factorizable operators can also contribute to the hadronization process. However, in the 
transverse/longitudinal and total e + e~ cross sections they are suppressed by additional fac- 
tors of l/Q and l/Q 2 , respectively, and thus usually negligible (Dasgupta and Webber, 
1997). 

The evolution of the fragmentation functions can be calculated by taking the logarithmic 
derivative of Eq. (5.2) with respect to the scale M 2 = Q 2 . This leads to the coupled 
homogeneous evolution equations (Georgi and Politzer, 1978; Baier and Fey, 1979; Altarelli 
et al, 1979) 



dD H/q (Q 2 ) _ a s {Q 2 
dlnQ 2 

dD H/g 



2tt 

+ P g ^ q e 

(Q 2 ) a s (Q 2 ) 



Pn 



dlnQ 2 



2tt 



„ - q & I) II ,,((/) 

D H/g (Q 2 ) 



Pq^g ® D H/q {Q< 



+ p, 



9^9 



d h/9 {Q 2 ) 



(5.4) 



where the time-like splitting functions Pj<-i are identical to the space-like splitting functions 
in LO (Altarelli and Parisi, 1977), but differ in NLO (Curci, Furmanski, and Petronzio, 
1980; Furmanski and Petronzio, 1980; Floratos, Kounnas, and Lacaze, 1981). 

The evolution equations can unfortunately not be solved analytically, and for light 
hadrons an initial distribution similar to the form Dh/%(x, Qq) = Nx a (l — x) 13 has to be 
assumed at a low starting scale Q . The free constants N, a, and (3 are then fitted to 
experimental data, usually from the process e + e~ — > j*(q) — > H(pn)X, where the initial 
state and the total center-of-mass energy Q 2 = q 2 are uniquely fixed. In the MS scheme, 
the normalized NLO semi-inclusive hadron cross section is given by (Baier and Fey, 1979; 
Altarelli et al, 1979) 



1 ^ Q2) -Y.A^) + ^ a ' m 



Ctot dx 7 , ; T± 



2tt 



(5.5) 



x {CI® D H/q (Q 2 ) + CI® D H/9 (Q 2 )}} , 



where 



1 , <*.(&) 



71 



(5.6) 
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Group 


Year 


Hadron 


Ql 

(GeV 2 ) 


Factor. 
Scheme 


N f 


Nm—A 
MS 

(MeV) 


BEP 


1979 




25.0 


LO 


3 


450/600 


AKL 


1983 


n ± ,K ± 


25.0 


LO 


3 


300/400 


GR+ 


1993/5 




900 


MS 


5 


269/319 


NW 


1994 




m% 


MS 


5 


344 


BKK 


1995 




2.00 


LO 


5 


190 










MS 




190 




1995/8 




2.00, 


LO 


5 


146 (fit) 






D *± B +,o 




MS 




317 (fit) 


CGRT 


1997 


D u , D u 


9 

mj 


MS 


5 


151 


KKP 


2000 


TT ± ,K ± , P 


2.00 


LO 


5 


121 (fit) 










MS 




299 (fit) 


Kretzer 


2000 


7T ± ,K ± ,h ± 


0.26 


LO 


5 


175 








0.40 


MS 




246 


BFGW 


2001 




2.00 


MS 


5 


300 



TABLE VIII. Parameterizations of various hadron fragmentation functions. 



is the total hadronic cross section, x = 2(pn • q)/Q 2 is the fraction of the center-of-mass 
energy transferred to the observed hadron H, and a = 4ira 2 / (3Q 2 ) is the total cross section 
for e + e~ — > At larger Q 2 , the couplings and propagators in Eqs. (5.5) and (5.6) are 

modified due to additional Z-boson exchange. The time-like transverse (T) and longitudinal 
(L) Wilson coefficients are given by (Altarelli et ai, 1979) 



Cg,T — C F 



3 3 1 

2 (1 " X) ~ 2(T^) _ 



1 + x 2 , 

+ 2 lnx 

1 — x 



+ (!+»») ~ X) 



+ 



2tt 2 9 



x 



+ 



5 





— Cf, 




= Cf 







{l-xf , . 2 . 2(1 
-ln[x 2 (l- x)] - 



x) 



X 



X 



X 



(5.7) 



The transverse cross section starts at O(a s ), while the longitudinal cross section starts only 
at 0(a s ). Thus in NLO the longitudinal 0(a 2 ) Wilson coefficients have to be included for 
longitudinal cross sections, but not for total cross sections (Rijken and van Neerven, 1996, 
1997). 

For a long time, only LO fragmentation functions for charged pions, kaons, and pro- 
tons were available, which had been fitted to e + e~ data from MARK I (Baier, Engels, and 
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Petersson [BEP], 1979) and had subsequently been updated using deep-inelastic data from 
EMC (Anselmino, Kroll, and Leader [AKL], 1983). Improved LO sets for pions, 77's, kaons, 
and D-mesons were obtained by comparing HERWIG Monte Carlo predictions in the next- 
to-leading logarithmic approximation with e + e~ data from TPC, PETRA, and later LEP 
(Greco and Rolli, 1993, 1995; Chiappetta et al, 1994; Greco, Rolli, and Vicini [GR+], 1995; 
Nason and Webber [NW], 1994; Cacciari et al. [CGRT], 1997). Recently several independent 
sets have been obtained by a full NLO evolution of a non-perturbative input at low starting 
scales and fits to precise e + e~ data from MARK II, TPC, PETRA, AMY, SLD, and LEP 
(Binnewies, Kniehl, and Kramer [BKK], 1995a, 1995b, 1996a, 1997, 1998a, 1998b; Kniehl, 
Kramer, and Pdtter [KKP], 2000a; Kretzer, 2000; Bourhis et al. [BFGW], 2001). These data 
make it possible to simultaneously determine a s with a precision that is competitive with the 
world average (Binnewies, Kniehl, and Kramer, 1995b; Kniehl, Kramer, and Potter, 2000b). 
The extracted value of a s could be in principle correlated with the size of power corrections 
(Nason and Webber, 1994), but the effect was found to be negligible (Kniehl, Kramer, and 
Potter, 2000b). According to the factorization theorem, fragmentation functions should be 
universal, i.e. process independent. The validity of this assumption has been verified by 
comparing inclusive hadron cross sections with e + e~ data not used for the fits at different 
scales Q 2 and with different hard scattering processes like photon-photon, photon-hadron, 
and hadron-hadron collisions (Binnewies, Kniehl, and Kramer, 1995b, 1996b, 1997, 1998a, 
1998b; Kniehl, Kramer, and Potter, 2001). The parameterizations of the various fragmen- 
tation functions discussed above are compared in Tab. VIII. From Fig. 26 one observes that 
the three NLO fits (KKP, Kretzer, and BFGW) are very similar for the fragmentation of 
light quarks and all quarks, except for the theoretically unreliable small-x region (x < 0.1) 
and the experimentally badly constrained large x-region (x > 0.8). For the fragmentation 
of c- and 6-quarks there are also differences at intermediate values of x. In the case of the 
6-quarks the discrepancy is due to the fact that two incompatible data sets from DELPHI 
(1991-93 and 1994 data) were used in the fits: Table IX shows that it is not possible to 
obtain good values of \ 2 P er degree of freedom, Xdf> simultaneously for both data samples. 
Otherwise all three NLO fits give satisfactory values of Xdf- 

While the masses of light quarks and gluons can always safely be neglected in the mass- 
less evolution approach, the masses of the heavier charm and bottom quarks h may 
become comparable to the physical scale Q. The even heavier top-quark decays before it 
hadronizes. Bjorken (1978) argued, that attaching a light antiquark q to a heavy quark h 
should decelerate the latter only slightly, so that it carries almost the same energy as the 
hadron H = (hq). As a consequence, a light quark fragmentation function proposed by Field 
and Feynman (1978), which was peaked at rather low x, was modified (Ali et al, 1980a, 
1980b). Later Peterson et al. (1983) expanded the energies in the fragmentation amplitude 
1 / (Eh — Eh + Eg) about the transverse particle masses with the result 

N 

DHIh{x) = x[l-l/x-e h /(l-x)r (5 ' 8) 

where oc m; \/rn\ is a free parameter and N is a normalization factor constrained by the 
momentum sum rule Eq. (5.3). This fragmentation function peaks at x — 1 — ^fEh and 
has a width ^fe\ l (Bodwin and Harris, 2001). Alternative forms have been proposed by 
Kartvelishvili, Likhoded, and Petrov (1978), Bowler (1981), Andersson et al. (1983), and 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



x 

FIG. 26. Dependence of the normalized inclusive hadron cross section at y/S = 91.2 GeV on the 
scaled hadron momentum x. Data from DELPHI (triangles for 1991-93 and circles for 1994 data) 
and SLD (squares) are compared with the NLO fragmentation functions of KKP (solid), Kretzer 
(dotted), and BFGW (dashed). Contributions from b, c, light, and all quarks (from bottom to 
top) are rescaled by factors of 1/5 relative to the nearest upper distribution (Kniehl, Kramer, and 
Potter, 2001). 
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Energy 


Flavor 


Experiment 


FF Set 


No. of 


(GeV) 






KKP 


Kretzer 


BFGW 


Points 


29 


uds 


TPC 


0.178* 


0.159 


0.167* 


7 




c 




0.876* 


0.911 


0.923* 


7 




b 




2.23* 


1.21 


1.14* 


7 


91.2 


all 


DELPHI 94 


1.28 


1.51* 


1.49 


12 






SLD 


1.32 


0.370 


0.421 


21 




uds 


DELPHI 91-3 


3.17* 


0.990* 


1.95 


13 






DELPHI 94 


0.201 


0.588* 


1.00* 


12 




c 


DELPHI 91-3 


0.473* 


0.388* 


0.401 


11 




b 


DELPHI 91-3 


28.9* 


0.887* 


1.03 


12 






DELPHI 94 


0.433 


9.14* 


8.74 


12 


189 


all 


OPAL 


0.568* 


0.250* 


0.414* 


11 



TABLE IX. Xdf va hies obtained in comparisons of the NLO fragmentation functions by KKP, 
Kretzer, and BFGW to e + e~ data with x > 0.1 at different center-of-mass energies. Data samples, 
which did not enter the respective fits, are marked by asterisks (Kniehl, Kramer, and Potter, 2001). 



Suzuki (1986). Using Heavy Quark Effective Theory, Braaten, Cheung, and Yuan (1993a) 
calculated the 0(a 2 s ) fragmentation function 

D Bc /b(x, m b + 2m c ) = N rX ^ ~ X \ [6 - 18(1 -2r)x 

[1 — (1 — r)x\° 

+ (21-74r + 68r 2 )x 2 - 2(l-r)(6-19r + 18r 2 )x 3 

+ 3(l-r) 2 (l-2r + 2r 2 )x 4 ] (5.9) 

for B c mesons, which depends on the mass ratio r = m c / (m b +m c ). Like the non-perturbative 
fragmentation functions discussed above, this result can also be used to parameterize the 
input at the starting scale Qo for heavy-light mesons with the normalization A" and r as 
free parameters (Braaten et al, 1995; Binnewies, Kniehl, and Kramer, 1998b). Eq. (5.9) 
peaks at large x (x ~ 0.8) and vanishes correctly like (1 — x) 2 as x — > 1, whereas a previous 
perturbative calculation by Collins and Spiller (1985) only vanishes like (1 — x) and fails to 
describe recent ALEPH data (Heister et a/., 2001). 

Matching a massless NLO calculation for e + e~ — > hX in the MS scheme to a NLO 
calculation with as a regulator for the collinear singularity, Mele and Nason (1991) 
calculated the functions 



D H/h (x,Q 2 ) = 5(l-x) + 



a s (Q 2 )C F 



l+x 2 / Q 2 

1 ln ~ 

1 — x \ mi 



2tt 

- 21n(l-z)-l)] + , 
D Blg (x, Q 2 ) = ^!^V + (1 - x f) In 
D H/q ^ h (x,Q 2 ) = 0(a 2 s ), (5.10) 
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which reproduce the massive cross section when convolved with the massless hard scattering 
cross section at the scale Q. For scales Q ^> m^, however, the logarithms in Eqs. (5.10) 
become large and have to be resummed. In this case, these equations can be interpreted 
as boundary conditions at a starting scale Q and can be evolved to the scale Q. Both 
approaches were found to give similar results in photoproduction (Cacciari and Greco, 1996) 
and e + e~ annihilation (Nason and Oleari, 2000), but they yielded a satisfactory description 
of the data only when an additional non-perturbative input Du/ i {x 1 Q'^) = Nx a {l — x) 13 
(Cacciari et al, 1997) or of the form in Eq. (5.8) (Cacciari and Greco, 1997; Nason and 
Oleari, 1999, 2000) was included. An alternative solution consists in adding Eqs. (5.10) to 
the massless cross section calculation with the interpretation of a change of factorization 
scheme (Kniehl, Kramer, and Spira, 1997; Binnewies, Kniehl, and Kramer, 1998a; Kramer, 
1999). The fitted values of the non-perturbative input depend, of course, on the treatment 
of the perturbative component. 

The fragmentation of quarks and gluons into bound states of heavy quarks and anti- 
quarks Q = (hh) is based on the general factorization analysis of Bodwin, Braaten, and 
Lepage (1995). It allows for the fragmentation functions of quarkonia Q to be factor- 
ized into calculable short-distance coefficients df(x, Q 2 ), which describe the production rate 
of a quark-antiquark pair with quantum numbers n = [a 2S+1 Lj] within a region of size 
1/m/i, and long-distance operator matrix elements (OMEs) (0 Q (n)), which contain the 
non-perturbative dynamics responsible for the formation of the bound state Q from the 
state n, 

D Q/l (x,Q 2 ) =J2dUx,Q 2 )(0 Q (n)}, ie{q,h,g}. (5.11) 

n 

The size of the OMEs may be estimated by how they scale with the relative velocity v of 
the quarks inside the quarkonium state (Lepage et al, 1992). In the case of the charmonium 
state J/^I the leading OME 

(O J 1 /q '( 3 S 1 )) = ^\RjM0)\ 2 (5-12) 

is related to the J/^> radial wave function at the origin Rj/q,(0) (Braaten, Fleming, and 

Yuan, 1996). The subleading OMEs (O^^S^), (O^fSo)), and (0 J 8 ^{ 3 Pj)) follow at the 
relative order 0(v 4 ). Of particular importance is the transverse gluon coefficient (Braaten 
and Yuan, 1994; Braaten and Lee, 2000; Ma, 1995) 

df^(x, Q 2 ) = ^1 U -x) + ^ (5.13) 

*[(a(-£ + ") + §-t+ 8 " 

+ 6 In 2 2) 5(1 - x ) + (\n^-- 1) P g ^ 9 (x) 
2 W1 ^ 6 ( ln(l — x) 



+ 6(2 - x + x 2 ) ln(l - x) - - 



x \ 1 — x 
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FIG. 27. Factorization of photon-hadron scattering into hadrons. 



since it starts already at 0(a s ), while the gluon coefficient for the leading OME (O-J ( 3 £i)) 
starts only at 0(al) (Braaten and Yuan, 1995). Since in the charmonium system v 2 ~ 
a s (m 2 ) ~ 0.3 (Quigg and Rosner, 1979), the (0$^( 3 Si)) contribution is no longer sup- 
pressed with respect to the (0(^ ( 3 Si)} contribution. The perturbative coefficient for the 
fragmentation of heavy quarks into the leading [l, 3 Si] state is of 0(a 2 ) (Braaten, Cheung, 
and Yuan, 1993b; Ma, 1994). 

For the P-wave Xc,j states, the coefficients of the leading OMEs (Of c " 7 ( 3 Pj)) with J = 
0, 1, 2 start at a 2 for both gluons (Braaten and Yuan, 1994; Ma, 1995) and quarks (Ma, 
1996a). The coefficients for the subleading OMEs (0£ c,J ( 3 Si)) start at 0(a s ) for gluons (see 
above) and 0(a 2 ) for quarks (Ma, 1996a), as do the quark (Braaten, Cheung, and Yuan, 
1993b; Ma, 1994) and gluon (Braaten and Yuan, 1993; Ma, 1994) coefficients for the leading 
OME in r] c production (O^^Sq)). 



B. Light hadrons 

Like jets, hadrons can be produced in direct, single-resolved, and double-resolved photon- 
photon scattering and in direct and resolved photon-hadron scattering. As can be seen in 
Fig. 27 for the photoproduction case, there is a new factorization in the final state for 
the transition of partons c into observed hadrons H, which is described by fragmentation 
functions D H / C . In the hadronic cross section 

/ dx a dx b — f a/A (x a , M 2 )f b/B (x b , M 2 ) 



&p 2 T <ly 



xD H/c (z,M 2 )^, (5.14) 



the fragmentation function is convolved at a factorization scale M c with the partonic cross 
section Eq. (4.10) and the initial state parton densities f a /A and fb/B, which were already 
encountered in the jet production case. 
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FIG. 28. Transverse momentum (left) and rapidity (right) dependence of the NLO inclusive 
charged hadron photoproduction cross section (Kniehl, Kramer, and Potter, 2001) compared to 
HI data (Adloff et ai, 1999a). 



The potential of inclusive hadron production for distinguishing the direct and resolved 
contributions was stressed already in early LO calculations for photon-photon (Brodsky et 
ai, 1979) and photon-hadron scattering (Fontannaz et ai, 1980). For light hadrons, the 
calculation of Born, virtual loop, and real emission processes proceeds as described in Sec. 

IV. B. There is, however, an important difference. Since the degeneracy of the final state is 
incomplete due to the observation of an identified hadron, the collinear singularity associated 
with the real emission of a massless parton 3 from the observed final parton line 2* 

\M F \ 2 ab ^ 123 = -- £ P^Ax)\M B \ 2 ab ^ + O(e ) (5.15) 

remains uncancelled and has to be factorized into the fragmentation function (see Sec. 

V. A). NLO calculations based on the formalism of Ellis et ai. (1980) have been performed 
for direct (Khalafi, Landshoff, and Stirling, 1983; Aurenche et ai, 1985a; Gordon, 1994), 
single-resolved (Aurenche et ai, 1984b, 1987; Gordon, 1994), and double- resolved (Aversa 
et ai, 1989) light hadron production, and their dependences on the renormalization and 
factorization scales were found to be considerably reduced with respect to the LO calcu- 
lations (Borzumati, Kniehl, and Kramer, 1993; Kniehl and Kramer, 1994; Gordon, 1994; 
Greco, Rolli, and Vicini, 1995). The NLO calculations were then used to demonstrate the 
universality of NLO fragmentation functions for charged hadrons, neutral pions, and neu- 
tral kaons (Binnewies, Kniehl, and Kramer, 1995b, 1996a, 1996b; Binnewies, 1997; Kniehl, 
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FIG. 29. Transverse momentum (left) and rapidity (right) dependence of the NLO inclusive neutral 
kaon photoproduction cross section (Kniehl, 1997) compared to HI data (Adloff et al, 1997). 



1997; Kniehl, Kramer, and Potter, 2001) by comparing them with photoproduction data 
from HI (Abt et al, 1994; Linsel, 1995; Adloff et al, 1997, 1999a) and ZEUS (Derrick et 
al, 1995c; Breitweg et al., 1998f) and photon-photon data from TASSO (Brandelik et al., 
1981), MARK II (Cords et al, 1993), L3 (Achard et al, 2001), and OPAL (Ackerstaff et 
al, 1999). 

In Fig. 28, NLO calculations of the inclusive charged hadron photoproduction cross 
section at HERA with AFG photon and CTEQ5M (Lai et al, 2000) proton densities and 
three different choices of the scale factor £ = /i/pr = M 7 /p T = M p /p T (Kniehl, Kramer, and 
Potter, 2001) are compared to HI data (Adloff et al, 1999a). The transverse momentum 
spectrum (left) agrees well with the data in shape and normalization. The rapidity spectrum 
(right) is sensitive to variations of the photon densities (AFG, GRV, or GS96), particularly 
for the lower p T cut, but the scale uncertainty is still large. 

A similar comparison is made in Fig. 29 for the photoproduction of neutral kaons. Due 
to the larger kaon mass, the massless NLO calculation describes the measured pt spectrum 
(left) only at larger px (pr > 1-5 GeV). The cut on p? for the rapidity spectrum (right, pr>2 
GeV) may still be too low for quantitative conclusions, but the potential to constrain f g / 7 
and D K o/ g is clearly visible. 

NLO charged hadron cross sections in photon-photon collisions (Kniehl, Kramer, and 
Potter, 2001) are shown in Fig. 30 and compared to OPAL data (Ackerstaff et al, 1999). 
Apart from the lowest pt points, which are below the starting scale of the fragmentation 
functions, the massless calculations describe the transverse momentum spectra (left) quite 
well. The scale variations in the rapidity spectra (right) are again significantly larger than 
the experimental errors and the sensitivity to the photon structure (not shown). 
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Pt [GeV] |y| 
FIG. 30. Transverse momentum (left) and rapidity (right) dependence of the NLO inclusive 
charged hadron cross section in photon-photon collisions (Kniehl, Kramer, and Potter, 2001), inte- 



grated over 77 invariant-mass intervals 10 < W < 30 GeV, 30 < W < 55 GeV, 55 < W < 125 GeV, 
and 10 < W < 125 GeV (from bottom to top in this order) and compared to OPAL data (Ackerstaff 
et al., 1999). 



C. Heavy hadrons 

If the transverse momentum px of a produced hadron is much larger than its mass, 
the formalism described in the last Section can be applied also to the production of heavy 
hadrons. The heavy quark is then considered to be one of the active flavors in the evolution 
of parton densities (see Sec. Ill) and fragmentation functions (see Sec. V.A), and the heavy 
quark mass is neglected in the hard scattering cross section. It is only kept as the starting 
scale in the evolution of the heavy quark parton densities and fragmentation functions. NLO 
calculations in the massless scheme have been performed for photoproduction (Cacciari and 
Greco, 1996, 1997; Kniehl et al, 1995; Kniehl, Kramer, and Spira, 1997; Binnewies, Kniehl, 
and Kramer, 1997, 1998; Kramer, 1999) and photon-photon collisions (Cacciari et al, 1996; 
Kramer and Spiesberger, 2001) with different assumptions on the perturbative and non- 
perturbative components of the fragmentation function (see Sec. V.A). They have been 
compared with D* photoproduction data from HI (Aid et al, 1996b) and ZEUS (Breitweg 
et al, 1997, 1999b; Derrick et al, 1997b) and photon-photon data from L3 (Acciarri et al, 
1999a) and OPAL (Abbiendi et al, 2000). 

If, on the other hand, the transverse momentum p T is of the same order as the mass of 
the produced hadron, the latter has to be kept in the partonic cross section 
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FIG. 31. Born diagrams for direct heavy-quark production. 
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FIG. 32. Born diagrams for single-resolved heavy-quark production. 

da B 1 1 / 4tts \ e 1 

~dT ~ 2^r(i-e) v(t-pi)(w-pi) 8^ 

x J^l^l 2 . 



S a SbC a Cb 



(5.16) 



In this case, the heavy flavors are not active in the parton densities and fragmentation 
functions, and only the Born diagrams in Figs. 31, 32, and Fig. 33 for direct, single-, and 
double-resolved photon-photon scattering contribute. The corresponding matrix elements 
are collected in Tab. X. Figs. 32 and 33 and the respective matrix elements apply also to 
direct and resolved photoproduction. In LO calculations (Fritzsch and Streng, 1978) one 
can of course set e = 0, but not if NLO corrections 



da v ^ 1 



4tts 



dt 2sF(l 

a. 



(t - pl){u - pi) -p 2 2 s 



8ns 



x 



r(i-e) 



S a S b C a C b 2n \ Q 2 / r(l-2e) 



\M 



V,F,I\2 



(5.17) 



are included. The heavy-flavor mass has to be kept in the denominators of the virtual loop 
integrals and in the integrations over the solid angle (Eq. (4.12)) and invariant mass 
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FIG. 33. Born diagrams for double-resolved heavy-quark production. 
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of the unresolved final state particles. The presence of a mass usually complicates the 
integrals that have to be computed, but the singularity structure is generally simpler than 
in the massless case, since the collinear singularities are regulated by the heavy-flavor mass. 
The method of analytic integration was originally developed for total and single-differential 
cross sections of heavy-flavor production in hadron collisions (Dawson, Ellis, and Nason, 
1988, 1989; Beenakker et al, 1989, 1991), but it was soon applied to photon-hadron (Ellis and 
Nason, 1989; Smith and van Neerven, 1992) and photon-photon collisions (Kiihn, Mirkes, 
and Steegborn, 1993; Drees et al, 1993; Contogouris, Kamal, and Merebashvili, 1995). Later 
different methods based on numerical cancellations of divergencies (Mangano, Nason, and 
Ridolfi, 1992) and phase-space slicing (Giele, Keller, and Laenen, 1996) were designed for 
more differential cross sections in hadron collisions. They were also subsequently extended to 
photon-hadron (Frixione et al, 1994a, 1994b, 1995; Frixione, Nason, and Ridolfi, 1995) and 
photon-photon collisions (Kramer and Laenen, 1996; Frixione, Kramer, and Laenen, 2000a, 
2000b). The massless limits of the massive NLO results have recently been compared with 
massless NLO calculations (Cacciari, Frixione, and Nason, 2001; Kramer and Spiesberger, 
2001). 

After convolution with various non-perturbative fragmentation functions, the massive 
NLO calculations were confronted with measurements of total charm cross sections, IP- 
meson px and y distributions, and DD azimuthal correlations by the Photon Emulsion 
(Adamovich et al, 1980, 1987), NA14/2 (Alvarez et al, 1992, 1993), E687 (Frabetti et al, 
1993, 1996; Moroni et al, 1994), E691 (Anjos et al, 1989, 1990), HI (Aid et al, 1996b), and 
ZEUS (Breitweg et al, 1997, 1999b, 2000e) photon-hadron experiments and by the JADE 
(Bartel et al, 1987), TASSO (Braunschweig et al, 1990), TPC/27 (Alston-Garnjost et al, 
1990), AMY (Takashimizu et al, 1996), TOPAZ (Enomoto et al, 1994a, 1994b), ALEPH 
(Buskulic et al, 1995), L3 (Acciarri et al, 1999b, 2001a, 2001b), and OPAL (Abbiendi et 
al, 2000, 2001a) photon-photon experiments. 

The massive NLO total charm cross section (Frixione et al, 1994a) shown in Fig. 34 
(top) describes the old fixed-target and the new HERA photoproduction data well within 
the scale uncertainty (Aid et al, 1996b). The HI measurement (Adloff et al, 1999b) of 
the total bottom cross section in Fig. 34 (bottom) is, however, severely underestimated by 
the massive NLO calculation, as are the corresponding ZEUS measurement (Breitweg et 
al, 2001) and the hadroproduction measurements from CDF (Abe et al, 1993, 1995) and 
DO (Abbott et al, 2000a, 2000b). It should be noted that the NLO predictions are rather 
sensitive to the photon (GRV) and proton (Martin, Roberts, and Stirling, 1993) parton 
densities, and to a lesser degree also to the fragmentation model and the masses of the 
charm (m c = 1.5 GeV) and bottom quarks (m b = 4.75 GeV), but these variations cannot 
account for the discrepancy in the bottom cross section (Blair et al, 2000). NLO massless 
calculations (Binnewies, Kniehl, and Kramer, 1997; Cacciari and Greco, 1997) cannot be 
applied to total cross sections due to the absence of a hard scale. They can, however, be 
applied to D* transverse momentum (top) and pseudorapidity (bottom) distributions as 
shown in Fig. 35. Massless distributions are usually larger than their massive counterparts, 
and they describe the experimental data better (Breitweg et al, 1999b), particularly in the 
large pt region, where the massive calculation only fits the data, if low values for m c = 1.2 
GeV, the renormalization scale /x = 0.5 • itlt, and the Peterson fragmentation function 
parameter e = 0.02 are chosen. The shapes of the measured pseudorapidity distributions 
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FIG. 34. Photon-proton center-of-mass energy dependence of the massive NLO total charm (top) 
and bottom (bottom) cross sections (Frixione et al, 1994a) compared to fixed target and HERA 
data (Aid et al, 1996b; Adloff et al, 1999b). 
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(Frixione et al., 1994a) and two different massless (Binnewies, Kniehl, and Kramer, 1997; Cacciari 
and Greco, 1997) NLO D* photoproduction cross sections compared to ZEUS data (Breitweg et 
al, 1999b). 
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are not well described by either the massive or the massless calculations. This fact has 
attracted considerable theoretical effort, from the consideration of "drag effects" due to the 
beam remnants (Norrbin and Sjostrand, 2000) to cc vector meson components of the photon 
wave function (Berezhnoy and Likhoded, 2001) and heavy-light quark recombination models 
inspired from heavy quarkonia (Berezhnoy, Kiselev, and Likhoded, 2000; Braaten, Jia, and 
Mehen, 2001), which may also have relevance for asymmetries in 5-physics. The cut on 
Pt > 2 GeV is probably too close to m c for the massless calculations to be reliable. In the 
remaining p? regions, different charm densities in the photon allow for a somewhat improved 
description of the data (Breitweg et al, 1999b). 

The measured total charm cross sections in photon-photon collisions in Fig. 36 (top) are 
again well described by massive NLO calculations (Drees et al, 1993), whereas the total 
bottom cross sections (bottom) are significantly underestimated. Even though the experi- 
mental and theoretical uncertainties are quite large, they cannot account for the discrepancy 
in the bottom case. As for the D* transverse momentum (top) and pseudorapidity (bottom) 
distributions shown in Fig. 37 they are well described by the massless calculation, despite 
the low transverse momenta studied {jpr > 2 GeV). Surprisingly, the massive calculation 
agrees with the data for larger pt {pt > 3 GeV), but not for small pt, where it should be 
more appropriate. Even low values of m c and /i cannot bring the massive calculation into 
agreement with the measured pseudorapidity distribution. 

With future e + e~ and ep colliders, it will not only be possible to extend the kine- 
matic ranges and production rates of charm and bottom quarks (Kramer and Laenen, 1996; 
Jankowski, Krawczyk, and Wing, 2001), but also to study the heavier top quark and its prop- 
erties (Kiihn, Mirkes, and Steegborn, 1993; Bernreuther, Ma, and McKellar, 1995; Choi and 
Hagiwara, 1995). 

D. Quarkonia 

As already briefly discussed in Sec. V.A, the factorization formalism of Non-Relativistic 
QCD (NRQCD) (Bodwin, Braaten, and Lepage, 1995) provides a rigorous theoretical frame- 
work for the description of heavy- quarkonium production and decay. This formalism im- 
plies a separation of short- distance coefficients, which can be calculated perturbatively as 
expansions in the strong-coupling constant a s , from long-distance operator matrix elements 
(OMEs), which must be extracted from experiment. The relative importance of the latter 
can be estimated by means of velocity scaling rules (Lepage et al, 1992), i.e. they are 
predicted to scale with a definite power of the heavy-quark velocity v in the limit w < 1. 
In this way, the theoretical predictions are organized as double expansions in a s and v. A 
crucial feature of this formalism is that it takes into account the complete structure of the 
hh Fock space, which is spanned by the states n = [a, 2S+1 Lj] with definite spin S, orbital 
angular momentum L, total angular momentum J, and color multiplicity a — 1, 8. In 
particular, NRQCD predicts the existence of color-octet processes in nature. This means 
that hh pairs are produced at short distances in color-octet states and subsequently evolve 
into physical (color-singlet) quarkonia by the nonperturbative emission of soft gluons. In 
the limit v — > 0, the traditional color-singlet model (CSM) (Berger and Jones, 1981; Baier 
and Riickl, 1981) is recovered. The greatest triumph of NRQCD was that it was able to 
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FIG. 37. Transverse momentum (top) and pseudorapidity (bottom) dependence of the massive 
(Frixione, Kramer, and Laenen, 2000a) and a massless (Binnewies, Kniehl, and Kramer, 1996b, 
1998a) NLO D* photon-photon cross section compared to OPAL data (Abbiendi et al, 2001a). 
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FIG. 38. Born diagrams for direct color-singlet J/* + 7 production. The charm quark lines in the 
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FIG. 39. Born diagrams for single-resolved color-singlet J/^+jet production. 



correctly describe the cross section of inclusive charmonium hadroproduction measured in 
pp collisions at the Fermilab Tevatron (Braaten and Fleming, 1995; Cho and Leibovich, 
1996a, 1996b), which had turned out to be more than one order of magnitude in excess of 
the theoretical prediction based on the CSM (Abe et al, 1992, 1997). An alternative model, 
which also predicts the existence of color-octet processes, is the duality or color evaporation 
model (CEM), where conservation of the color quantum number is ignored in the hard in- 
teraction and assumed to be recovered only through the exchange of non-perturbative soft 
gluons with the underlying event (Fritzsch, 1977; Halzen, 1977; Gluck, Owens, and Reya, 
1978). The short-distance color-singlet and color-octet cross sections are related to those 
of open heavy-quark production by statistical factors and to those of the physically ob- 
served mesons by constant universal factors. While this model can also account for the 
Tevatron data (Amundson et al, 1997), it lacks the status of NRQCD as an effective field 
theory. Similar objections may be raised to a model, which allows to flip spin and color 
of the heavy-quark pair through a perturbative gluon exchange with an isotropic classical 
color field generated by bremsstrahlung of the initial state partons (Hoyer and Peigne, 1999; 
Hoyer, Marchal, and Peigne, 2000). This "hard comover scattering" model can explain not 
only the observed ratio of Xc,i an d Xc,2 hadroproduction cross sections, but also the observed 
non-polarization of hadroproduced J ftp and ip' mesons. In both of these cases the CSM and 
NRQCD have some difficulty in explaining the data. Quarkonium production in pp and ep 
collisions has recently been reviewed by Kramer (2001). 

In order to convincingly establish the phenomenological significance of the color-octet 



72 



g ORnr* 



=h(T(r(p g g of(f(r*==(x>= J/XJ/ snf(n 




j/¥ 



uTCP g 



FIG. 40. Born diagrams for double-resolved color-singlet J/^+jet production. 
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TABLE XL LO matrix elements squared |.M B | 2 for direct, single-resolved, and double-resolved 
color-singlet J/* production. 



processes, it is indispensable to identify them in other kinds of high-energy experiments 
as well. In photon-photon and photon-hadron collisions, quarkonia can be produced by 
direct, single-resolved, and double- resolved processes. Which of the various photon, quark, 
and gluon initiated processes contribute to the production of a given quarkonium state, 
depends on the leading OMEs in the NRQCD velocity expansion and on the additional 
particles in the final state. In inelastic processes, the quarkonium state is produced with 
non-zero transverse energy and a balancing jet or photon. In this case the inclusive hadronic 
production cross section is given by Eq. (5.14) with Dh/ c {z, M%) = 5(1 — z). px and y now 
relate to the transverse momentum and rapidity of the observed quarkonium state, and 
the partonic cross section is given by Eq. (5.16), where the mass p\ is identified with the 
quarkonium mass M 2 (= Am^ in the static limit). 

Since its discovery in 1974, the J/^f charmonium has remained in the center of interest 
due to its relatively large production cross section and easy experimental identification in 
the leptonic decay modes. The diagrams contributing to the production of its leading color- 
singlet state [I, 3 Si] are shown in Figs. 38, 39, and 40. In photon-photon scattering, they 
relate to direct, single- resolved, and double-resolved processes. In photoproduction, the last 
two apply to direct and resolved processes, while in hadroproduction only the diagrams in 
Fig. 40 contribute. The corresponding hard scattering matrix elements are related by simple 
transformations of couplings and color factors and are listed in Tab. XI. The diagrams in 
Figs. 38, 39, and 40 contribute also to the production of one of the leading color-octet states, 
[8, 3 Si], when the final photon in Fig. 38 is replaced with a gluon and the color factors are 
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changed appropriately. At the same order v 4 relative to (0 1 J/ *( 3 1 S 1 )), the OMEs (O^Y^o)) 
and {O'l^ ( 3 Pj)} with J = 0, 1, 2 contribute in photon- and gluon-quark scattering with final 
state quarks and in quark-antiquark annihilation with final state gluons (Ma, McKellar, and 
Paranavitane, 1998; Japaridze and Tkabladze, 1998; Godbole, Indumathi, and Kramer, 
2001; Klasen, Kniehl, Mihaila, and Steinhauser, 2001a). 

In the photon and gluon fusion processes described above, the two internal quark prop- 
agators are off-shell by p\, so that the differential cross section falls like f/px- On the other 
hand, when pt ^> M, the quarkonium mass M can be considered small and the cross sec- 
tion should scale like any other single-particle inclusive cross section as l/p?- The dominant 
quarkonium production mechanism at sufficiently large pp must therefore be via fragmenta- 
tion of a final state parton (see Sec. V.A) in the partonic processes shown in Figs. 7, 8, and 
9 (Braaten and Yuan, 1993). 

In Fig. 41, the LO color-singlet and the sum of color-singlet and color-octet contributions 
are compared with preliminary data on J /iff production in photon-photon collisions from 
DELPHI (Todorova-Nova, 2001; Chapkin, 2002). The agreement of the NRQCD prediction 
with the data is excellent, although the statistical and theoretical uncertainties from the 
renormalization and factorization scales, from m c = 1.5 GeV, and the numerical values 
of the NRQCD operator matrix elements (Braaten, Kniehl, and Lee, 2000) are still large 
(Klasen, Kniehl, Mihaila, and Steinhauser, 2001b). Single- resolved processes are far more 
important than double-resolved and direct processes (see also Godbole, Indumathi, and 
Kramer, 2001). The total color-singlet contribution is insufficient to describe the data, and 
the size of the cross section is sensitive to the color-octet operator values extracted at the 
Tevatron. The shapes of color-singlet and color-octet contributions can, however, only be 
distinguished at larger p?. 

At the scale of the charm mass, the strong coupling a s is still relatively large (~ 0.3) and 
NLO corrections can become important. The real corrections to J /iff production in direct 
photon-photon collisions, which become dominant at large pt, have recently been calculated 
with an invariant mass cut-off for the associated dijet system. Fig. 42 shows that the LO 
cross section is enhanced by about an order of magnitude at large Pt, where the missing 
virtual corrections are no longer important and the uncertainty from the exact value of the 
cut-off is reduced (Klasen, Kniehl, Mihaila, and Steinhauser, 2001a). In a complete NLO 
calculation, there would of course be no dependence on the cut-off. 

A full NLO calculation is available for inelastic color-singlet J /iff production in direct 
photon-hadron collisions (Kramer et ai, 1995; Kramer, 1996). The direct process dominates 
at the relatively large values of z (z — pp-pj/y/pp-p-y — Ej/^/Ej > 0.3) currently observed 
at HERA. Recent HI data (Kriiger, 2000) show evidence also for resolved photoproduction, 
which is closely connected to hadroproduction and contributes mostly at low z < 0.3. This 
region is dominated by color-octet processes in NRQCD (Beneke, Kramer, and Vanttinen, 
1998; Kniehl and Kramer, 1999), but also in the CEM (Eboli, Gregores, and Halzen, 1999). 
As can be seen in Fig. 43, the direct LO color-singlet prediction can almost account for 
the data. Within the relatively large uncertainty of the theoretical parameters (shaded 
band), the direct color-octet processes improve the agreement at intermediate z, but show a 
stronger rise than the data at large z (Cacciari and Kramer, 1996; Ko, Lee, and Song, 1996; 
Amundson, Fleming, and Maksymyk, 1997). However, this region is close to the kinematic 
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FIG. 41. LO squared transverse momentum distribution of J/^Sl mesons produced in photon- 
photon collisions at LEP2. The total NRQCD result describes the preliminary data from DELPHI 
(Todorova-Nova, 2001; Chapkin, 2002) well, while the CSM is an order of magnitude too low. 
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FIG. 42. Transverse momentum spectrum for J/*S> production in direct photon-photon collisions 
at LEP2. At large px, the NLO real emission corrections enhance the LO cross section by an order 
of magnitude (Klasen, Kniehl, Mihaila, and Steinhauser, 2001a). 

limit and is influenced by the non-perturbative soft transition of the color-octet state into 
the physical J meson, so that LO predictions are not reliable. NLO corrections for color- 
octet photoproduction have unfortunately only been calculated for the total cross section 
(Maltoni, Mangano, and Petrelli, 1998), but higher-order terms have been resummed using 
non-relativistic shape functions (Beneke, Rothstein, and Wise, 1997; Beneke, Schuler, and 
Wolf, 2000). Relativistic corrections to the color-singlet state have been evaluated, but there 
is unfortunately no agreement on their effects in the literature (Greub et al., 1993; Khan 
and Hoodboy, 1996; Ma, McKellar, and Paranavitane, 2000). Higher twist effects due to 
multiple interactions with the proton or photon remnants should be suppressed as a power 
of A 2 / (Q 2 + pj.) with Q = m c , m c v, or m c v 2 (Beneke, Kramer, and Vanttinen, 1998) and 
not as a power of A 2 /(4m 2 ) (Ma, 1997), since only the first result leads to the intuitively 
expected A 2 /p^ behavior at large pr- On the other hand, the discrepancy between the 
NRQCD prediction and the HERA data can clearly be reduced, if higher order effects of 
multi-gluon emission are approximately taken into account in the extraction of the operator 
matrix elements at the Tevatron (Kniehl and Kramer, 1999). 

While the shape of the J/^> energy distribution in Fig. 43 is already well described by 
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FIG. 43. Direct and resolved contributions to the color-singlet and color-octet J/^f energy distri- 
bution at HERA (Kramer, 2001). 



the LO color-singlet calculation, the shape of the LO transverse momentum distribution 
in Fig. 44, which falls like a 2 s m^./p T , differs markedly from the HERA data. The data 
clearly require the NLO color-singlet contribution, which is dominated by t-channel gluon 
exchange and falls only like a^m 2 /p T (Kramer et al., 1995; Kramer, 1996). A similar scaling 
is observed by contributions from charm quark fragmentation in direct photoproduction, 
which exist at 0(al). They are suppressed at lowpr, but exceed the color-singlet production 
via photon-gluon fusion at px > 10 GeV according to LO (Godbole, Roy, and Sridhar, 
1996) and NLO calculations (Kniehl and Kramer, 1997a, 1997b). Gluon fragmentation in 
resolved photoproduction exceeds the fusion process already at px > 5 GeV due to a large 
enhancement by NLO corrections (Kniehl and Kramer, 1997a, 1997b). 

Further information about the underlying production mechanism (CSM, CEM, or 
NRQCD) may be gained from the energy spectrum of J/\l/'s produced in polarized photon- 
hadron collisions (Chao et al, 2000; Japaridze, Nowak, and Tkabladze, 2000) and of J/\l/'s 
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FIG. 44. J/ty transverse momentum distribution in LO and NLO direct color-singlet photopro- 
duction (Kramer, 2001). 



produced in association with a photon, where the color-octet contribution dominates already 
at intermediate z (Kim and Reya, 1993; Cacciari, Greco, and Kramer, 1997; Mehen, 1997). 
Another very sensitive observable is the energy spectrum of the polar angle in polarized 
J/ty decays (Beneke, Kramer, and Vanttinen, 1998). At the upgraded HERA collider one 
may furthermore hope to observe and study the photoproduction of other quarkonia like 
the S'-wave states ip' (Morii, Roy, and Sudoh, 2000; Brugnera, 2001) and r\ c (Chao, Hao, 
and Yuan, 1999) or the P-wave states \ c (Ma, 1996b) and h c (Fleming and Mehen, 1998). 
Bottomonia have a smaller relative squared velocity v 2 ~ 0.1 than charmonia with v 2 ~ 0.3 
(Quigg and Rosner, 1979) and should therefore be better described by the NRQCD velocity 
expansion. Unfortunately, the photoproduction cross section for T is more than two orders 
of magnitude smaller than the one for J/ty due to the smaller electric charge and larger 
mass of the bottom quark. 

The large conceptual and numerical uncertainties in quarkonium physics clearly call 
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for an increased theoretical effort, in particular for a full NLO calculation of color-singlet 
production at hadron colliders and of color-octet production at photon and hadron colliders. 
In open flavor production, the NLO corrections should be calculated for the photoproduction 
of two light or heavy hadrons and of hadrons in association with jets, which would help 
to disentangle the flavor (e.g. charm) content of photons and protons and to resolve the 
question, why massive calculations cannot describe the observed bottom cross sections. 



79 



VI. PROMPT PHOTON PRODUCTION 



The many facets, that photons exhibit in the initial state of hard scattering processes, 
can be rediscovered in their production in the final state: If produced inclusively, prompt 
photons are related to hadrons by important fragmentation contributions; if produced in 
isolation, they exhibit more perturbative aspects due to their pointlike coupling to quarks; 
if produced in association with jets, they become powerful tools for constraining the initial 
parton densities. All of these aspects will be discussed in this Section. 



A. Fragmentation 



Since photons couple to charged quarks, they can be produced directly in hard scattering 
interactions. If the photon is radiated from a massless final-state quark, it exhibits a collinear 
1/e- singularity proportional to the P^ q (x) splitting function and the underlying LO matrix 
element. This 1/e-pole has to be absorbed into a non-perturbative fragmentation function, 
just like the collinear singularity of an initial-state photon had to be absorbed into a non- 
perturbative parton density. The photon fragmentation functions D y /i are thus intimately 
related to the parton densities in the photon /j/ 7 . 

The evolution equations of these fragmentation functions are of a form similar to Eqs. 
(3.5) with the role of initial and final particles reversed. They are given by (Koller, Walsh, 
and Zerwas, 1979) 
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for a single quark flavor q. The time-like splitting functions Pj^i are identical to the space- 
like splitting functions in LO (Altarelli and Parisi, 1977), but differ in NLO (Curci, Furman- 
ski, and Petronzio, 1980; Furmanski and Petronzio, 1980; Floratos, Kounnas, and Lacaze, 
1981). Like the evolution equations Eqs. (3.5) for photon densities, Eqs. (6.1) contain inho- 
mogeneous terms: The quark-photon splitting function can be obtained in LO from P g ^ q 
by the transformation P^ q = e 2 JCpPg^q- The gluon-photon splitting function enters only 
in NLO and is given by (Aurenche et al, 1993) 
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Information about the photon fragmentation function can be obtained from the process 
e + e~ — > 7X. However, this process also has a large contribution from the direct coupling of 
photons to quarks, which dominates over the fragmentation contribution, particularly if the 
photon is experimentally isolated (see Sec. VLB). The direct contribution has been used 
extensively to extract the electroweak couplings of up- and down-type quarks to the Z-boson 
by the ALEPH (Decamp et al, 1991; Buskulic et al, 1993), DELPHI (Abreu et al, 1995), 
L3 (Adriani et al., 1992, 1993), and OPAL collaborations (Akrawy et al., 1990; Alexander 
et al, 1991; Acton et al, 1993). 

In the MS scheme the NLO inclusive photon cross section is given by (Altarelli et al, 
1979) 



a dx q 9 { ' 2tt q 7 2tt 

x [c T q ® D l/q {Q 2 ) + Cj ® D l/g {Q 2 )] } , (6.3) 

where a = Ana 2 N c / (3Q 2 ) is N c times the pointlike cross section for e + e~ — > , 
the factor of two comes from D^/ q (Q 2 ) = D 1 / q (Q 2 ), and the time-like Wilson coefficients 
for the transverse and longitudinal cross sections are given by Eqs. (5.7). The photonic 
Wilson coefficients are = C^JCf. C^ t contains a term oc ln[x 2 (l — x)] with a loga- 
rithmic singularity at x — 1 and a stronger singularity at x — 0. Gliick, Reya, and Vogt 
(1993) defined a modified DIS 7 scheme for the photon fragmentation function by absorb- 
ing the transverse photonic Wilson coefficient CT uT into the quark fragmentation function 

D^iQ 2 ) = D ^f q {Q 2 ) + a /( 2 ^)e 2 q Cl T . This affects again the NLO photon splitting func- 
tions 

P^=P^_ q -e 2 q P q ^ q ®Cl T , 

= P™ "2E e ? P^ 9 ® CIt ■ (6-4) 

In LO of a, the third evolution equation in Eqs. (6.1) can be directly integrated with the 
result D 7 / 7 (x, Q 2 ) — 5(1 — x). Furthermore in LO of a s only the evolution equation of the 
quark-photon fragmentation function 



dD l/q (x,Q 2 ) _ <x_ 
dlnQ 2 ~ 2tt 



P^ q {x) (6.5) 



survives, which can also be integrated with the result (Gehrmann-de Ridder, Gehrmann, 
and Glover, 1997) 

a Q 2 

D l/q (x, Q 2 ) = —P^ q (x) ln gf + D l/q {x, Q ). (6.6) 

The first term in Eq. (6.6) is the perturbatively calculable pointlike solution, while the 
second term is a hadronic boundary condition which has to be fitted to experimental data, 
e.g. on e + e~ — > 7+jet. Following this approach the ALEPH collaboration obtained the LO 
result (Buskulic et al, 1996) 
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D l/q (x, Ql) = — (-P^ q (x) ln(l - xf - 13.26) 



(6.7) 



at the starting scale Qo = 0.14 GeV. For the fragmentation function of quarks to photons 
with virtuality P 2 , the purely perturbative LO result is (Qiu and Zhang, 2001; Braaten and 
Lee, 2002) 




(6.8) 



The NLO evolution equation for the quark-photon fragmentation function can also be inte- 
grated, if one neglects the higher-order contributions from the gluon-photon fragmentation 
function, the mixing between different quark flavors, and the running of a s (Q 2 ). The result 
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Ql 



+ D y/q (x, Ql) 



is then exact at the fixed order 0(aa s ), but does not take into account the usual resummation 
of powers of ln(Q 2 /Ql). ALEPH have fitted their e + e~ — > 7+jet data in NLO, taking 
a s = 0.124 in order to reproduce the observed total rate of e + e~ — > hadrons, and obtain 
(Buskulic et a/., 1996) 



D l/q (x, Ql) = (-P^ q (x) ln(l - x) 2 + 20.8(1 - x) 



a 

- 11.07) 



(6.10) 



at Qo = 0.64 GeV. 

The fully resummed solution of the evolution equations consists of pointlike and hadronic 
contributions D^ fi (Q 2 ) = D^ fi (Q 2 ) + D^jf(Q 2 ). The full LO (Duke and Owens, 1982) and 
NLO pointlike solutions (Aurenche et al, 1993; Bourhis, Fontannaz and Guillet, 1998) 
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pi 



a 



7 /i(Q ) - 27f 



47T 



a s (Q 2 ) 



di + bi + 0(a s 



(6.11) 



can only be calculated analytically in moment space. The hadronic input can unfortu- 
nately not be determined from inclusive photon production in e + e~ annihilation, since 
the experimental data (Buskulic et al, 1996; Ackerstaff et al, 1998) are very limited and 
furthermore dominated by the pointlike quark-photon fragmentation function. Therefore 
all existing parameterizations (Aurenche et al. [ACFGP] 1993; Gliick, Reya, and Vogt 
[GRV], 1993; Bourhis, Fontannaz, and Guillet [BFG], 1998) use VMD to model the photon 
fragmentation at low scales. The different input parameters are summarized in Tab. XII. 
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FIG. 45. Parameterizations of the photon fragmentation function at Q 2 = M§ compared to OPAL 
data (Ackerstaff et al, 1998). 
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Group 




Set 


Ql 
(GeV 2 ) 


Factor 
Scheme 


VMD 
Model 


Nt 

IV f 


(MeV) 


ACFGP 


1993 


NLO 


2.00 


MS 


coherent 


4 


230 


GRV 


1993 


LO 


0.25 


LO 


incoherent 


5 


200 






NLO 


0.30 


DIS 7 






200 


BFG 


1998 


NLO 


0.50 


MS 


coherent 


4 


230 



TABLE XII. Parameterizations of the photon fragmentation functions. The coherence of the VMD 
model is determined by the coefficients e u and in Eq. (3.16). ACFGP and GRV relate the photon 
to the pion, while BFG fit p production data from ALEPH and HRS. 



Heavy quarks are included above their production thresholds with boundary conditions 
D l/h (x,m 2 h ) = D^ rh {x,m 2 h ) = 0. 

Fig. 45 demonstrates that the most recent data from OPAL (Ackerstaff et al, 1998) can 
already be described by the pointlike fragmentation functions in LO (long-dashed curve) 
and NLO (dotted curve). The data can thus not discriminate between the different LO (full 
curve) and NLO (short-dashed and dot-dashed curves) assumptions for the hadronic input. 

B. Isolation 

Photons produced via fragmentation usually lie inside hadronic jets, while directly pro- 
duced photons tend to be isolated from the final state hadrons. The theoretical uncertainty 
coming from the non-perturbative fragmentation function can therefore be reduced if the 
photon is isolated in phase space. At the same time the experimental uncertainty coming 
from photonic decays of 7r°, 77, and u mesons is considerably reduced. Photon isolation can 
be achieved by limiting the (transverse) hadronic energy inside a cone of size R around 
the photon to 

e It) < £(t)E(t),t (6.12) 

This is illustrated in Fig. 46. For small photon pseudorapidities 77, the cone size R = 
5/ cosh 77 is approximately equal to the half angle of the cone 5. Since in NLO the direct and 
fragmentation processes are intimately linked at the factorization scale Q, it is mandatory 
that the isolation criterion does not interfere with the cancellation of soft and collinear 
singularities. For certain differential cross sections like da/dx^ in e + e~ annihilation, where 
x 1 = 2E 1 /Q 2 is the fractional photon energy, spurious infrared singularities remain at the 
point x 7 = 1/(1 + e) (Berger, Guo, and Qiu, 1996). The cross section da/dx^ is therefore 
only defined as a distribution, which has to be integrated with a test function over finite bin 
widths. So for physical observables these integrable singularities disappear. Nevertheless, 
the cross section contains logarithms of the type ln[l — x 7 (l + e)], which can become large 
at x 7 — 1/(1 + e) (Aurenche et al, 1997; Catani, Fontannaz, and Pilon, 1998). For physical 
observables it is therefore important to integrate over large enough bins in x 1 . The situation 
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FIG. 46. Illustration of an isolation cone containing a parton c that fragments into a photon 7 plus 
hadronic energy Ef Tag . In addition, a gluon enters the cone and fragments giving hadronic energy 
-Eparton (Berger, Guo, and Qiu, 1996). 

is very similar to the infrared sensitivity of the NLO dijet cross section at large values of the 
initial state photon energy fraction i* s . Recently an improved photon isolation criterion 



has been proposed, where 5 < 5q and 5q is now the isolation cone (Frixione, 1998a). This 
procedure allows the fragmentation contribution to vanish in an infrared safe way. 

C. Inclusive photons 

The photoproduction of inclusive photons proceeds through direct and resolved initial 
state photons and through direct and fragmentation production of the final state photon. 
Thus four types of partonic subprocesses contribute to 7/1 — > 7X, as shown in Fig. 47, and 
eight types to 77 — > 7X. The hadronic cross section is given by Eq. (5.14), where the hadron 
fragmentation function is replaced by the photon fragmentation function D 7 / c . For direct 
photon production the fragmentation function is, of course, simply 5(1 — z), and the cross 
section coincides with Eq. (4.14). 

The partonic LO diagrams and cross sections can be obtained by crossing the diagrams in 
Figs. 7, 8, and 9 and the matrix elements in Tab. IV, and the NLO corrections are calculated 




(6.13) 
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FIG. 47. Factorization of lepton-hadron scattering into photons. 




as outlined in Sec. V.B. The only additional subtlety occurs when a final state quark splits 
into a quark and a photon. The corresponding collinear singularity 
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(6.14) 



is absorbed into the photon fragmentation function 

D l/q (x, M 2 ) = D l/q (x) + [T^ q (M 2 ) ® D lh ] (x) (6.15) 

at the factorization scale Mj. 

NLO corrections to resolved-direct and resolved-fragmentation processes have been cal- 
culated in the context of prompt photon production in hadronic collisions by different groups 
(see Sec. VII. C), while those for the direct-fragmentation process have been calculated in the 
context of inclusive hadron photoproduction (see Sec. V.B). The NLO corrections, which 
are specific only for the photoproduction of prompt photons, i.e. the corrections to the 
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FIG. 48. The quark box contribution to direct photon photoproduction. 



direct-direct process, have been calculated for inclusive photons by Duke and Owens (1982), 
Aurenche et al (1984a), Bawa, Krawczyk, and Stirling (1991), and Gordon and Storrow 
(1994). Full NLO calculations for photon-hadron scattering have been performed by Gor- 
don and Vogelsang (1995), using a subtraction term in order to account also for photon 
isolation, and by Fontannaz, Guillet, and Heinrich (2001a) in a completely differential form. 
Gordon and Storrow (1996) have furthermore calculated NLO corrections to prompt photon 
production in photon-photon collisions. 

The direct process 77 — > 7X can only proceed through a quark box diagram with X = 7 
and is of 0(a A ). The process 75 — > 7X with X = g proceeds through the quark box shown 
in Fig. 48 and similar diagrams with reversed and crossed fermion flow and is of 0(a 2 a 2 ). It 
thus represents a NNLO correction to the tree-level subprocesses which contain at least one 
quark line and are of 0(a 2 ). Since the box contribution is by itself gauge invariant and can 
be numerically large, it can and has to be included together with the NLO corrections of the 
tree-level processes (Aurenche et al, 1984a; Krawczyk and Zembrzuski, 2001a; Fontannaz, 
Guillet, and Heinrich, 2001a). 

In Fig. 49 the only existing complete NLO calculation for isolated photons, which also 
includes the NNLO box contribution (Fontannaz, Guillet, and Heinrich, 2001a), is compared 
to recent data from ZEUS (Breitweg et al, 2000d). The transverse momentum (pj) distri- 
bution (top) agrees quite well with the data, but the experimental errors are still too large to 
distinguish between the AFG and GRV parton densities in the photon. On the other hand, 
the NLO pseudorapidity (r? 7 ) distributions with AFG photon densities (bottom), integrated 
over j? T e [5; 10] GeV, show some discrepancies with the data: The calculations underesti- 
mate the data for negative rapidities and small fractional photon energies y = E 7 /E e and 
overestimate them for positive rapidities and large photon energies. Similar results are ob- 
tained with the calculations of Gordon (1998) and Krawczyk and Zembrzuski (2001a), which 
are partly based on LO cross sections. The movement of the maximum towards the backward 
direction is caused by an increase of the photon energy participating in the hard scattering 
process. A possible explanation for the discrepancy in the forward direction might be that 
there is more hadronic activity in this experimental region, which is more strongly affected 
by the isolation cut than the theoretical parton level simulation. At a future THERA col- 
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FIG. 49. Transverse momentum (top) and pseudorapidity (bottom) spectra of isolated photons in 
NLO compared to ZEUS data (Fontannaz, Guillet, and Heinrich, 2001a). 
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lider, where 250 GeV electrons would be collided with 920 GeV protons, the cross section in 
the ZEUS kinematic region would be about three times larger (Krawczyk and Zembrzuski, 
2001b). 

D. Photons and jets 

Prompt photon production in association with an observed jet has been calculated by 
Gordon [LG] (1998) and by Krawczyk and Zembrzuski [KZ] (2001a). However, both of these 
calculations use LO expressions for the resolved-fragmentation contribution (the latter also 
for the direct-fragmentation and resolved-direct contributions). Numerical results of these 
calculations, using GRV and GS parton densities in the photon, are compared with ZEUS 
data (Lee et al, 2000) and PYTHIA (Sjostrand, 1994) and HERWIG (Abbiendi et al, 1992) 
Monte Carlo predictions in Fig. 50. The shape of the transverse energy distribution (top), 
integrated over rp G [—0.7; 0.9], is described by all theoretical predictions, but the magnitude 
is not quite reproduced by the Monte Carlo predictions. The same observation can be made 
in the pseudorapidity distribution (bottom) integrated over Ej, e [5; 10] GeV, but here the 
higher order calculations with GS photon densities also fall below the data in the backward 
region. The main difference between the two higher order calculations comes from the box 
contribution, which is included only in the calculation of Krawczyk and Zembrzuski (2001a). 

Fixed-target and hadron collider experiments indicate that prompt photon production 
may be influenced by intrinsic transverse momenta of the scattering initial partons (see 
Sec. VII. C). ZEUS have analyzed isolated prompt photon photoproduction in association 
with a jet and observed an excess over the LO QCD expectation in the distributions of 
perpendicular momentum p±, longitudinal momentum py, and azimuthal angle A0 of the 
photon relative to the balancing jet (Chekanov et al, 2001). They have then attributed 
these discrepancies to an effective (k T ) of 1.69 ±0.18to;2o GeV, which includes effects coming 
from the initial-state parton showering as modeled within PYTHIA (Sjostrand, 1994). This 
value of (k T ) seems to be consistent with determinations in hadron collisions at different 
energies, which are, however, obtained using a variety of methods. Within the LO Monte 
Carlo comparison, the ZEUS data seem to support the trend that the effective (kr) in 
the proton rises with the available hadronic energy. However, the data have recently been 
confronted with a calculation for photon plus jet photoproduction, that includes all NLO 
corrections and the NNLO box diagram (Fontannaz, Guillet, and Heinrich, 2001b). The 
comparison is complicated by the fact that equal cuts on the transverse energies of the 
photon and recoiling jet have been used in the experiment (Et > 5 GeV), which has to be 
relaxed in the NLO calculation and induces some sensitivity on the size of the mismatch 
l-^T,min ~~ ^T,minl < 0-5 GeV. The sensitivity within this range is, however, smaller than 
the experimental uncertainty. Within errors, the full NLO calculation can describe the 
ZEUS data without any intrinsic (kx) effects. The production of direct photons and jets in 
direct, single-, and double-resolved photon-photon collisions has also recently been calculated 
in NLO (Fontannaz, Guillet, and Heinrich, 2002), and good agreement with preliminary 
OPAL data (Abbiendi et a/., 2001d) was found. In the future, NLO calculations should 
be performed for the production of two photons and photons in association with light and 
heavy hadrons in order to obtain more information about the flavor content of photons and 
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FIG. 50. Transverse energy (top) and pseudorapidity (bottom) spectra of isolated photons pro- 
duced in association with a jet of Et > 5 GeV, —1.5 < r/ < 1.8, and cone size R = 1 compared to 
ZEUS data (Lee, 2000). 
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VII. RELATED TOPICS 



In this Section, the topics discussed so far will be extended in several directions. First 
the photon virtuality, which was assumed to be negligible so far, will be explicitly taken 
into account. Second, the initial photon will be allowed to have a well-defined polarization, 
and third the production of prompt photons in hadron collisions will be discussed, which is 
closely related to photoproduction by crossing an initial state photon into the final state. 



A. Virtual photons 



If the momentum transfer in a lepton scattering process q = pi — py is small, but not 
completely negligible, the corresponding cross section can be calculated in the Equivalent 
Photon approximation using a differential form of the transverse photon flux (Budnev et ai, 
1974; Kessler, 1975) 
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(7.1) 



The Weizsacker- Williams approximation (Eq. (2.1)) is recovered by integrating Eq. (7.1) 
over the photon virtuality P 2 = —q 2 . As P 2 becomes larger, the longitudinal photon flux 
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(7.2) 



also has to be taken into account. Bremsstrahlung spectra for transverse photons with 
different virtualities are shown in Fig. 51 (top) together with the ratio of longitudinal and 
transverse photon fluxes (bottom), which is essentially independent of the photon virtuality. 

Like real photons, photons with a non-zero virtuality P 2 can be probed for their hadronic 
structure in deep-inelastic electron-photon scattering. As P 2 approaches the virtuality Q 2 of 
the probing photon, the structure of the virtual photon reduces to the perturbative splitting 
of the photon into a collinear quark- ant iquark pair with probability (Uematsu and Walsh, 
1981, 1982; Rossi 1984) 



f^(x,Q 2 ,P 2 
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Q 2 
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(7.3) 



This expression includes contributions from transverse and longitudinal photons. It differs 
from the real photon case f g / y (x, Q 2 ) = e 2 aC^/ (2tt) with C 7 given in Eq. (3.11), since 
the collinear singularity is now mass regulated by the photon virtuality P 2 . In the limit 
P 2 <C Q 2 , the structure of the virtual photon should, of course, coincide with that of the 
real photon. 

A simple phenomenological ansatz to suppress smoothly the quark and gluon content in 
the virtual photon at high P 2 is (Drees and Godbole [DG], 1994) 
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Bremsstrahlung Spectra for Virtual Photons 
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FIG. 51. Bremsstrahlung spectra for transverse photons with different virtualities (top) and the 
ratio of longitudinal and transverse photon fluxes (bottom) . The latter is essentially independent 
of the photon virtuality. 
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where the typical hadronic scale P c lies between the QCD scale A and the proton mass. 
Since gluons are radiated from off-shell quarks, they must be more suppressed than quarks 
(Borzumati and Schuler, 1993). A similar ansatz has also been used by Aurenche et al. 
(1994b). It can be applied to any parameterization for the parton densities in the real 
photon, which are recovered in the limit P 2 <C Q 2 . At large P 2 , the evolved parton densities 
are globally suppressed. This means, however, that the perturbative result of Eq. (7.3) is 
not reproduced, since the logarithmic suppression factors are independent of x. 

A refined version of the above prescription has been formulated by Schuler and Sjostrand 
[SaS] (1995, 1996), who sum the real hadronic contributions from the three dominant vector 
mesons p, u, and after weighting them with factors f](P 2 ) = (l + P 2 / m p,o;,^) _2 and evolving 
them from max(P 2 , Ql). The pointlike contribution is evolved with a continuous suppression 
factor (1 + P 2 /k 2 )~ 2 , so that it also starts effectively from max(P 2 , Q 2 ,). In the low- and 
high-P 2 limits, the real photon densities and the logarithmically enhanced terms of the 
virtual pointlike box are correctly reproduced, but not the x-dependent terms of the latter. 

Gliick, Reya, and Stratmann [GRSt] (1995) interpret Eq. (7.3) at P 2 = Q 2 and 
f g /j(x, P 2 , P 2 ) = as NLO pointlike boundary conditions. In LO these conditions dif- 
fer since here also f q / 1 {x, P 2 , P 2 ) = 0. The pointlike boundary conditions are then added to 
the hadronic VMD input for real photons fffi(x, Qq), adopted from Gliick, Reya, and Vogt 
[GRV] (1992b), with an interpolating factor r?(P 2 ) = (1 + P 2 /m 2 )" 2 , so that 

f ih (x, max(P 2 , Ql), P 2 ) = V (P 2 )f-ff(x, max(P 2 , Q 2 )) 

+ [1 - v(P 2 )W*(x, max(P 2 , Ql), P 2 ). (7.5) 

The functional form of i](P 2 ) can be derived from a dispersion integral in the photon virtu- 
ally, which is dominated by the poles of the vector mesons, particularly the p-meson, in the 
region of small P 2 . In the limit P 2 Q 2 , the purely hadronic real photon input is recovered 
in LO and NLO, as required by the DIS 7 scheme. As P 2 approaches Q 2 , the unevolved 
initial parton densities are globally suppressed and the virtual pointlike result Eq. (7.3) is 
recovered. This is, however, inconsistent with the original DIS 7 scheme, since the NLO 
pointlike quark boundary condition differs now from zero, but is evolved using the same 
massless DIS 7 splitting functions (Eqs. (3.13)) in the inhomogeneous evolution equations 
(Eqs. (3.5)) from the new starting scale max(P 2 , Q 2 ,)- A parameterization of these virtual 
photon parton densities, which is valid in the region P 2 < Q 2 /5, has been performed only 
in LO (Gliick, Reya, and Stratmann, 1996). 

The inconsistencies with the DIS 7 scheme can be avoided if the pointlike term in Eq. 
(7.5) is omitted (Gliick, Reya, and Schienbein [GRSc], 1999b, 2001). However, in this case 
the perturbative splitting Eq. (7.3) is not reproduced. Consequently, this virtual parame- 
terization is only applicable in the region of small P 2 (P 2 < Q 2 /10), where the virtuality of 
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Parameterizations of the Partem Densities in the Virtual Photon 
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FIG. 52. NLO parameterizations of the up-quark (top), clown-quark (center), and gluon (bottom) 
densities in a photon with virtuality P 2 = 0.5 GeV 2 (left) and 2.5 GeV 2 (right) at Q 2 = 25 GeV 2 . 
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the pointlike term is of minor importance, and the large-P 2 region should be calculated in 
fixed order. A parameterization has also been given only in LO, although the NLO results 
are very similar. 

The existing NLO parton densities in the transverse virtual photon are compared in 
Fig. 52 at Q 2 = 25 GeV 2 and P 2 = 0.5 GeV 2 (left) and 2.5 GeV 2 (right) for the up- 
quark (top), down-quark (center), and gluon (bottom) densities. The suppression of the 
real photon GRV densities in the virtual photon is clearly visible. For the quarks (gluons) 
it is weakest (strongest) for the global DG factors, for which P 2 = 0.3 GeV 2 has been 
used. At P 2 = Q 2 /50 the ^-dependence is only weakly modified, while at P 2 = Q 2 /10 
the logarithmic singularity at x = in the perturbative box is significantly weaker. The 
different parameterizations vary largely at small P 2 , particularly for the gluon, but much 
less at larger P 2 as expected. 

Information about the parton densities in the virtual photon can be gained from the total 
hadronic cross section in double-tagged e + e~ scattering in the Bjorken region P 2 <C Q 2 . 
If the probing and target photons are both soft, their transverse (T) and longitudinal (L) 
polarizations contribute with equal weights to the effective virtual photon structure function 

F? s ~ F T \ + F T \ + i£ T + ift ~ Fl + \Fl (7.6) 

In the second equation, the relations F^ h ~ F^ T = Fl and F^ L ~ have been used, which 
apply strictly speaking only to the virtual box contribution with F2 = J2 q ^xe 2 -^e 2 N c Ax(l — 
x). Like F2 , Fl could be calculated by inserting longitudinal virtual photon densities into Eq. 
(3.12) and employing the longitudinal Wilson coefficients C Qj l = 8x/3 and C 9i l = 2x(l —x). 
In practice, however, only a pointlike parameterization for the longitudinal virtual photon 
is currently available (Chyla and Tasevsky, 2001), and it turns out to be very close to the 
unevolved longitudinal box result. 

In Fig. 53, F e g measurements by the PLUTO (top) (Berger et al, 1984b) and L3 (bottom) 
(Acciarri et al, 2000) collaborations are compared with the virtual photon parameteriza- 
tions discussed above. In order to obtain F^ s from F% , the separately shown longitudinal 
virtual box contribution 3/2F^ h has been added to all parameterizations. For the GRSc 
parameterization, which is evolved with three massless flavors, the sum of direct and re- 
solved heavy charm quark contributions is also shown separately. While the PLUTO and 
L3 data tend to be higher than all parameterizations and the purely perturbative virtual box 
prediction, the theoretical curves all describe the measurements within the errors, so that 
present data do not yet convincingly establish the importance of the renormalization group 
improved treatment of parton densities in the virtual photon or the need for an intrinsic 
hadronic component (Gliick, Reya, and Schienbein, 2001). 

The virtual photon fluxes and parton densities discussed above make it possible to extend 
the calculations of real photoproduction processes to the region P 2 > 0, as long as P 2 is sig- 
nificantly smaller than the hard scattering scale Q 2 (typically the squared transverse energy 
E\ of the observed jet, hadron, or prompt photon). Therefore virtual photoproduction will 
also receive contributions from direct and resolved photon processes. They can be calculated 
by taking P 2 into account only in the photon flux and parton densities, but neglecting it in 
the hard scattering. When P 2 approaches Q 2 , however, the resolved contribution vanishes 
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Parameterizations of the Effective Virtual Photon Structure Function 



PLUTO (Q 2 = 5 GeV 2 , P 2 = 0.35 GeV 2 ) 



Full Box 
GRSc (NLO) 
GRSt (NLO) 
SaS ID (LO) 
SaS 2D (LO) 




L3 (Q 2 = 120 GeV 2 , P 2 = 3.7 GeV 2 ) 



Full Box 
GRSc (NLO) 
GRSt (NLO) 
SaS ID = 2D (LO) 




FIG. 53. Parameterizations of the effective virtual photon structure function compared to data 
from PLUTO (top) (Berger et al, 1984b) and L3 (bottom) (Acciarri et al, 2000). 
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logarithmically, the longitudinal photon cross section has to be taken into account, and the 
DIS cross section has to be recovered. 

As already discussed in Sec. Ill, direct and resolved real photoproduction processes are 
related in NLO by the appearance of a collinear 1/e singularity in the splitting of the 
photon into quarks with fractional charge e q (Eq. (3.2)). This singularity is absorbed into 
the renormalized quark densities in the real photon (Eq. (3.3)) and cancels its logarithmic 
factorization scale dependence. A similar mechanism occurs for virtual photoproduction, 
where the singularity is now mass regulated by P 2 and Eq. (3.2) has to be replaced with 



This defines an MS factorization scheme for virtual photons, which leaves identical finite 
terms in the hard cross section as for real photons (Klasen, Kramer, and Potter, 1998; 
Potter, 1999c). Since there is no 1/e singularity, the subtraction is not mandatory to obtain 
finite cross sections, but only to resum the logarithm which can become large for P 2 <C Q 2 
and x — > 1. 

In order to obtain more information about the structure of the virtual photon it is useful 
to consider not only inclusive electron-photon scattering, but also hard scattering processes 
like jet photoproduction. LO studies indicate a significant contribution from resolved virtual 
photons at HERA, if P 2 <C E\ (Gliick, Reya, and Stratmann, 1996; Gliick, Reya, and 
Schienbein, 2001). Virtual jet photoproduction has also been calculated in NLO by adapting 
an existing NLO calculation for resolved real photons (Klasen, Kleinwort, and Kramer, 
1999) and using the MS subtraction for the direct virtual photon singularity discussed above 
(Klasen, Kramer, and Potter, 1998). It has then been extended to the DIS region of large P 2 
and applied to jet production at HERA (Kramer and Potter, 1998, 1999) and LEP (Potter, 
1999a, 1999b). The resolved contribution turns out to be particularly important, if the jets 
are produced in the forward region. The ZEUS (HI) collaborations have measured jets with 
rj > 2.6 (rj E [1.735; 2.794]) and E T > 5 GeV (E T > 3.5 GeV), respectively (Breitweg et al, 
1999c; Adloff et al, 1999c). In Fig. 54 the ZEUS data are compared to NLO predictions. 
While the DIS direct calculation (a) significantly underestimates the ZEUS data, they can 
be described if resolved virtual photon contributions from the SaS ID fit are added (b) 
(Kramer and Potter, 1999). 

In principle, the production of heavy quarks with mass is also sensitive to the struc- 
ture of the virtual photon. In this case resolved photon contributions are applicable as long 
as P 2 <C m 2 h . However a LO study shows that they are numerically small at the relatively 
large values of x currently probed at HERA (Gliick, Reya, and Stratmann, 1996) and LEP 
energies (Gliick, Reya, and Schienbein, 2001). 

Finally, a LO study for virtual photoproduction of prompt photons has shown that the 
gluon content in the virtual photon may be tested if the photons are produced in the forward 
direction at transverse energies E T ^ 5 GeV (Krawczyk and Zembrzuski, 1998). 



T q ^ r (x,P 2 ) = 6 <n 6(l-x) 
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P g ^(x) + N c e 2 q 



(7.7) 
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FIG. 54. NLO predictions for forward dijet production with virtual photons compared to ZEUS 
1995 data. In (a) only the DIS direct contribution is included, while in (b) direct and resolved 
virtual photon contributions have been added consistently (Kramer and Potter, 1999). 
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B. Polarized photons 



Polarized photons can be produced by bremsstrahlung off lepton beams with circular 
polarization |Aj| < 1/2. The corresponding spectrum Af^f ms (x) = /^ ems ' + (x) - f^™' ( x ) 
is given by (Philipsen, 1992; De Florian and Frixione, 1999) 

'1 - (1 -xf 1 2m 2 x 2 ' 



A A f brems 

dP 2 ^^)~ 27r 



X 



p2 



p4 



(7.8) 



for transverse virtual photons and by 
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(7.9) 



for almost real photons. As in the unpolarized spectrum (Eq. (2.1)) a non-logarithmic 
term is present, which is, however, not singular for x — > 0. At x ~ 1 the photons are 
completely polarized parallel to the incoming lepton helicity, but at x ~ 0, where most of 
the bremsstrahlung photons are produced, they are completely unpolarized (see dot-dashed 
curve in Fig. 55). The lepton polarization is lost in the Lorentz transformation from the 
Breit frame of the lepton-photon vertex to the center-of-mass frame of the photon-target 
vertex. 

The circularly polarized beamstrahlung spectrum has only been derived in an integral 
form (Schroeder, 1990) 
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where the Airy-function Ai(t>) falls off exponentially at large v and where u = 
[5x] 2 / 3 /[4 v / 3T(l -x)] 2 / 3 . The polarization Af^ m (x)/ f^ m (x) depends only weakly on the 
beamstrahlung parameter T and is very similar to the bremsstrahlung polarization (Berge, 
Klasen, and Umeda, 2001). 

While the photon polarization at an e + e~ collider is thus rather limited, a photon collider 
offers the additional possibility to control the helicity of the laser photons \P C \ < 1. The 
outgoing photons have a polarized spectrum (Ginzburg et al, 1984) 
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FIG. 55. Degree of circular polarization for laser-backscattered photons. The dot-dashed curve 
applies also to brems- and beamstrahlung photons, but then it extends out to x = 1 and 

\/ef{x)/f 1 / e {x) = 1. 



If only the electrons are polarized (P c = 0), the x-dependence of the polarization 
^f^fe V ( x )l f X ^fe T ( x ) i s si m ilar to the brems- and beamstrahlung cases (dot-dashed curve in 
Fig. 55). However if P c = ±1, then the helicity of the backscattered photons is opposite to 
that of the laser photons at x = x m & x (full, dashed, and dotted curves in Fig. 55). Therefore 
the choice 2\ e P c = —1 guarantees not only good monochromaticity (see Fig. 1), but also 
a high degree of polarization of the produced photons. By switching the signs of A e and 
P c simultaneously, one can switch the helicity of the outgoing photons without spoiling the 
monochromaticity of the photon spectrum. 

The polarized parton densities A/j/ 7 (Q 2 ) of the photon obey the same perturbative evo- 
lution Eqs. (3.5) as their unpolarized counterparts, except that polarized splitting functions 
APj^i have to be used. These have recently been calculated in NLO by Mertig and van 
Neerven (1996) and by Vogelsang (1996a, 1996b). We review only the LO results (Altarelli 
and Parisi, 1977) 



C F 



q^q 



l + x 2 3 r/i . 
—) + + 2 6{1 - X) 
(*), 



+ 0{a s 



101 



AP, 



g^q 



» 



'l-(l-x) 



X 



x 2 -(l-x) 2 } +0{a s ), 
1 + x 4 ) 



v z (l-x)J 



X 



+ 



6 ° 



3 1 



8(l-x) + 0(a s ), 



(7.12) 



where the Dirac matrix 75 has been assumed to anticommute with the Dirac matrices 7 M . 
Different results are obtained in other 75-schemes as mandated in dimensionally regularized 
NLO calculations. Some of them require finite renormalizations to arrive at the correct 
final answer. The polarized photon-quark splitting function can be obtained in LO by 
the transformation AP g ^_ 7 = 2Nce 2 AP q <_ g , and the polarized NLO photon-gluon splitting 
function is given by (Stratmann and Vogelsang, 1996) 
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As in the unpolarized case, the pointlike solution of the polarized evolution equations 
dominates at large x and Q 2 . It takes the same functional form as Eq. (3.9), where Aai and 
Abi replace and hi and have been explicitly calculated in moment space in LO (Irving and 
Newland, 1980; Hassan and Pilling, 1981; Xu, 1984) and NLO (Stratmann and Vogelsang, 
1996). The position of the low-a; singularities in Aa« and Abi differs from the unpolarized 
case, since the polarized splitting functions are now involved. 

By combining the polarized parton distributions functions Af q / 1 (Q 2 ) and Af g / 1 (Q 2 ) with 
the appropriate Wilson coefficients (Ratcliffe, 1983; Bodwin and Qiu, 1990) 
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and (Stratmann and Vogelsang, 1996) 
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one obtains the NLO photon structure function in the MS scheme 
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The ln(l — x)-term in the Wilson coefficient of the polarized photon causes similar stabil- 
ity problems in g\ as in F^ at large x. Thus it should again be absorbed into the quark 
distributions, which affects also the polarized pointlike NLO splitting functions. The trans- 
formation is given by Eqs. (3.13) with all quantities replaced by their polarized counterparts. 
This DISa 7 factorization scheme leads then also to perturbatively stable results and to purely 
hadronic boundary conditions in LO and NLO. 

For unpolarized photons, the hadronic input can be inferred from pionic parton densities. 
These are, however, not known in the polarized case. Thus Stratmann and Vogelsang (1996) 
were forced to make the two extreme assumptions that either Af^(Q 2 j ) = fffi(Ql) or 
Afffi(QQ = 0, as had already been suggested in the LO analysis of Gliick, Stratmann, and 
Vogelsang [GStV] (1994). Current conservation implies that the first moment of g± vanishes 
(Bass, 1992; Narison, Shore, and Veneziano, 1993), which can be realized by demanding 
Af^' ra_1 = 0. This is trivial to fulfill in the second (minimal) scenario, but requires 
additional assumptions on the low-x behavior of the first (maximal) scenario. In a previous 
LO analysis, Gliick and Vogelsang (1992) considered these two scenarios only for the gluon, 
while they fixed Afffi(Ql) = fffi(Ql)- The original DIS 7 scheme, which was constructed 
to avoid negative values of i 7 ^ at large x, unfortunately causes a violation of the positivity 
constraint \gj\ < F± = (i^ — Fj)/{2x). For polarized photons it is thus preferable to 
absorb the photonic Wilson coefficient of F?, Ci )7 = C 7 — AN c x{l — x), into the pointlike 
splitting functions Eqs. (3.13) (Gliick, Reya, and Sieg [GRSi], 2001a). The maximal polarized 
boundary conditions in the DISa 7 scheme then have to be identified with the hadronic 
input in the DIS 7i i scheme, and the minimal polarized boundary conditions are given by 
Affj^(Ql)e 2 q a/ {2it)(C\ n — C 7 ). All the polarized photon densities discussed above are based 
on the unpolarized GRV parameterization (Gliick, Reya, and Vogt [GRV], 1992b), which 
imposes an incoherent VMD input at a low starting scale and is evolved in moment space 
with up to five massless flavors (see Sec. III.C). The different polarized parton densities in 
the photon are compared in Fig. 56. 

Polarized photon densities can also be extended to non-zero photon virtualities, either 
by multiplying the boundary conditions with a suppression factor i](P 2 ) = (1 + P 2 /m 2 )~ 2 
(Gliick, Reya, and Sieg, 2001b) to guarantee a smooth transition to real photons, or by 
imposing the polarized virtual box result (Sasaki and Uematsu, 1999, 2000, 2001) 

A/{£(s, Q 2 , P 2 ) = 3^(2* _ i) (l„ ^_ _ 2 ) (7.17) 

at Q 2 = P 2 . As in the unpolarized case, the virtual box can be suppressed at small P 2 by 
a factor 1 — i](P 2 ). 

The polarized cross sections for single-, two-, and three-jet photoproduction can be ob- 
tained from Eqs. (4.14), (4.21), and (4.27) by replacing the unpolarized photon spectra, 
parton densities, and partonic cross sections with their polarized counterparts Af 7 /i(x), 
^fih(x,Q 2 ), and 

dAa B 1 1 2g 2 



dt 2s Ms S a S b C a C b 



(|A*?+| 2 -|A<?_| 2 ), (7.18) 



where parity invariance guarantees that |.M+ + | 2 = |.M B _| 2 and |.M+_| 2 = |A^^ + | 2 . Simi- 
larly the polarized cross section for the photoproduction of light hadrons can be obtained 
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Parameterizations of the Parton Densities in the Polarized Photon 
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FIG. 56. Parameterizations of the polarized up-quark (top), down-quark (center), and gluon (bot- 
tom) densities in the photon in LO (left) and NLO (right) at Q 2 = 25 GeV 2 . 
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from Eq. (5.14). The massless Born diagrams are the same as those in Figs. 7, 8, and 9. 
The purely partonic (double-resolved) diagrams were evaluated by Babcock, Monsay, and 
Sivers (1979) using the projection operators 



with an arbitrary four-vector n (n 2 = 0) for quarks and gluons with momentum p and 
definite helicity h. From their results, the direct and single-resolved matrix elements can 
easily be obtained by appropriate changes of couplings and color factors. They can then be 
applied to jet production in a polarized collider mode of HERA (Stratmann and Vogelsang, 
1997) or to inclusive hadron photoproduction in fixed-target collisions (Airapetian et ai, 
2000; Contogouris, Grispos, and Veropoulos, 2000). Polarized NLO calculations have been 
performed for direct photoproduction of inclusive hadrons (De Florian and Vogelsang, 1998) 
and also for resolved (De Florian et al, 1999) and direct (De Florian and Frixione, 1999) 
photoproduction of one or two jets. If large luminosities could be accumulated in a polarized 
HERA or eRHIC collider, these processes could be used to determine the polarized parton 
densities in the photon and proton, in particular A/ fl / 7 and Af g / p . In this way the infor- 
mation obtained from polarized structure functions in DIS would be supplemented. Fig. 57 
shows the pseudorapidity dependence of the single-jet asymmetry Aa/a at HERA energies. 
The asymmetry is clearly sensitive to the polarized photon (top) and proton (bottom) den- 
sities. In addition, as a ratio of cross sections it is very stable with respect to higher order 
corrections and variations of the factorization and renormalization scales. 

Polarized heavy quark photoproduction proceeds through the Born diagrams shown in 
Figs. 31, 32, and 33, which applied already to the unpolarized production. Gliick and Reya 
(1988) have calculated the photon-gluon fusion squared matrix element, which is also appli- 
cable to the photon-photon process after changing couplings and color factors. The resolved 
squared matrix elements have been evaluated by Contogouris, Kamal, and Papadopoulos 
(1990) and by Karliner and Robinett (1994), but they do not contribute significantly at the 
energies available at HERA and a possible GSI collider with E e = 5 GeV and E p = 50 GeV 
(Stratmann and Vogelsang, 1997). Due to this reduced sensitivity to the unknown polarized 
structure of the photon, heavy quark photoproduction may be a useful tool to constrain the 
polarized gluon density in the proton, despite the fact that fragmentation effects have to be 
taken into account (Frixione and Ridolfi, 1996). As expected the theoretical uncertainties 
are considerably reduced if NLO corrections for the direct (Contogouris, Kamal, and Mere- 
bashvili, 1995; Jikia and Tkabladze, 1996) and single-resolved (Bojak and Stratmann, 1998, 
1999; Contogouris, Grispos, and Merebashvili, 2000a, 2000b) processes are included. For 
the double-resolved process the NLO corrections are not yet available. 

NLO corrections are also unavailable for polarized quarkonium production in the color- 
singlet and in the color-octet channels. However LO studies suggest that polarization pro- 



u(p, h)u(p, h) 



v(p, h)v(p, h) 




(7.19) 
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FIG. 57. Pseudorapidity dependence of the single-jet asymmetry in polarized photoproduction for 
different polarized photon (top) and proton (bottom) densities (De Florian and Frixione, 1999). 
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FIG. 58. Factorization of hadron-hadron scattering into photons. 

vides an even more stringent test of NRQCD factorization than unpolarized quarkonium pro- 
duction, since the color-octet asymmetries are significantly different from the color-singlet 
asymmetries and since the uncertainties from the quarkonium operator values cancel in the 
asymmetries to a large extent. On the other hand, quarkonium production can also be used 
to constrain the polarized gluon density in the proton, provided that the color-octet operator 
values can be fixed in a different place (Chao et al, 2000; Japaridze, Nowak, and Tkabladze, 
2000; Morii and Sudoh, 2000). 

C. Prompt photons in hadron collisions 

The production of prompt photons in hadron collisions is intimately related to photopro- 
duction by crossing the initial state photon into the final state. The contributing hadronic 
scattering processes are shown schematically in Fig. 58, where the photon is either produced 
directly in the hard partonic subprocess (left) or via fragmentation of a final state parton 
(right). The inclusive cross section for the direct process Eq. (4.14) is a convolution of the 
parton densities in the initial hadrons, which can be pions, protons, antiprotons, or nuclei, 
with the partonic cross section. The fragmentation cross section Eq. (5.14) contains an ad- 
ditional convolution with the photon fragmentation function D 7 / c . The LO parton diagrams 
and matrix elements for the fragmentation process are the same as those in Fig. 9 and Tab. 
IV, while those for the direct process are obtained by crossing the photon leg in Fig. 8 and 
Tab. IV. In these processes the photon is balanced by an outgoing jet or hadron. Double 
prompt-photon production can be calculated by taking into account the direct diagrams and 
matrix elements in Fig. 7 and Tab. IV and a second photon fragmentation function for the 
double-resolved processes discussed above. 

The calculation of the NLO corrections proceeds along the lines discussed in Sec. V.B. 
For the inclusive direct process, they have been evaluated by Aurenche et al. (1984c, 1988, 
1990) for unpolarized and by Gordon and Vogelsang (1993, 1994a, 1994b) also for polarized 
prompt photon production. The production of a direct photon in association with an ob- 
served jet has been calculated in NLO by Baer, Ohnemus, and Owens (1990a, 1990b), by 
Gordon (1997a, 1997b), and by Frixione (1998b, 1999). The last two calculations have again 
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FIG. 59. Prompt photon data from fixed-target and ISR experiments as a function of the scaled 
transverse momentum xt and normalized to NLO QCD predictions with ABFOW proton parton 
densities and BFG photon fragmentation functions (Aurenche et ai, 1999). 
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included beam polarization. Polarized and unpolarized direct photon production with an 
additional final state charm quark has been calculated in the massless scheme by Berger and 
Gordon (1996, 1998) and by Bailey, Berger, and Gordon (1996). The direct NLO corrections 
to double prompt-photon production have been evaluated by Aurenche et al. (1985b), Bai- 
ley, Ohnemus, and Owens (1992), Coriano and Gordon (with polarization) (1996a, 1996b), 
and Binoth et al. (2000, 2001). NLO corrections to the fragmentation process have been 
evaluated by Aversa et al. (1989) and were applied to prompt photon production by Au- 
renche et al. (1993). Fragmentation contributes typically less than 20 % in fixed-target 
collisions, but can become dominant at collider energies. An isolation cut can then help to 
significantly reduce the theoretical and experimental uncertainties from the fragmentation 
process (Berger and Qiu, 1990, 1991). 

At large values of transverse momentum px, the LO QCD Compton process qg — > q^ 
dominates over the competing qq — > #7 annihilation process. This makes prompt photon 
production particularly sensitive to the gluon density in the proton, which enters only in NLO 
in DIS. Prompt photon production has indeed been used to constrain the unpolarized gluon 
density in the proton at large x (Vogelsang and Vogt, 1995; Martin et al, 2000) and could also 
be used to constrain the polarized gluon density in the proton (Gordon and Vogelsang, 1996; 
Gordon, 1997c; Chang, Coriano, and Gordon, 1998). Unfortunately, these determinations 
suffer from poorly known fragmentation contributions, particularly of the gluon and at small 
x (see Sec. VI. A), and potentially large logarithms in the isolation criterion (see Sec. VLB). 
Even worse, the data from fixed-target and ISR experiments cannot be described consistently 
by NLO calculations (Aurenche et al, 1999). As Fig. 59 demonstrates, particularly the 
data from the Fermilab E706 experiment lie a factor of two or more above NLO predictions 
(Apanasevich et al., 1998). They can only be described, if an intrinsic transverse momentum 
(k T ) ~ 1.3 GeV of the incoming partons is taken into account. This effect can either be 
generated by a simple Gaussian distribution or, as shown in Fig. 60, by considering multiple 
soft-gluon radiation and simultaneously resumming the large logarithms at small values of 
transverse momentum and at the partonic threshold xt = 2pr/VS ~ 1 (Catani, Mangano, 
Nason, Oleari, and Vogelsang, 1999; Laenen, Sterman, and Vogelsang, 2000). At collider 
energies the disagreement with NLO QCD is less severe and exists only at low transverse 
momenta, but it still leaves room for speculation on intrinsic (k^) effects (Abe et al, 1994; 
Abbott et al, 2000c). Due to the various uncertainties discussed above, prompt photon 
data have currently been dropped from global determinations of the gluon density at large 
x (Martin et al, 2001; Lai et al., 2000). Instead high-i? T jet data are used, which in turn 
suffer from a large factorization scheme dependence (Klasen and Kramer, 1996c; Anandam 
and Soper, 2000). 

Like real photon production, the production of photons with virtuality Q is also dom- 
inated by the QCD Compton process at large transverse momenta px > Q/2 and thus 
sensitive to the proton's gluon distribution (Berger, Gordon, and Klasen, 1998). However, 
the photon virtuality acts as a mass regulator and eliminates the fragmentation contribu- 
tions. Isolation is also not necessary, since the virtual photon can be identified by its muonic 
decay, and the presence of the scale Q may reduce the importance of large logarithms in pt- 
In Fig. 61 we show that virtual photon production in fixed-target collisions is indeed very 
sensitive to the gluon density (Berger and Klasen, 2000). In a similar way the polarized 
gluon distribution could be determined (Berger, Gordon, and Klasen, 2000a, 2000b). 
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VIII. SUMMARY 



From the fixed-target experiments in the 1980s to modern e + e~ and ep colliders, hard 
photoproduction has been a fertile research ground experimentally and theoretically. Precise 
measurements and calculations at next-to-leading order of perturbative QCD have been 
performed for a large variety of processes, and they have served a multitude of purposes. 
For one thing, the factorization theorem for hard scattering processes has been thoroughly 
tested. The photon energy spectra have been improved to meet the precision mandated by 
current and future lepton colliders. The parton densities of the photon, proton, and pion, and 
in particular the gluon densities, are now better constrained, which reduces the theoretical 
uncertainty at hadron and photon colliders. Ambiguities associated with jet algorithms and 
kinematic regions of large hadronization corrections have been identified and eliminated. 
Perturbative and non-perturbative aspects of parton fragmentation into individual hadrons 
and photons and of quarkonium formation have been disentangled, and the universality of 
fitted fragmentation functions and quarkonium operator expectation values has been tested. 

In spite of all these successes, research into hard photoproduction is far from complete: 
There is still large room for improvements, not only in even more precise next-to-next-to- 
leading order calculations, but also in determinations of the gluon density and spin structure 
of photons and protons, of the transition from real to virtual photons, and of photon in- 
teractions with new hypothetical particles as predicted e.g. by Supersymmetry. Exciting 
possibilities may soon open up with a new linear e + e~ collider, which may even include a 
dedicated photon-photon experiment with backscattered laser beams. Such an experiment 
would certainly initiate a whole new era of hard photoproduction. 
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